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Abstract 

We investigate the macroscopic dynamics of sets of an arbitrary finite number of 
weakly amplitude- modulated pulses in a multidimensional lattice of particles. The 
latter are assumed to exhibit scalar displacement under pairwise, arbitrary-range, 
nonlinear interaction potentials and are embedded in a nonlinear background field. 
By an appropriate multiscale ansatz, we derive formally the explicit evolution equa- 
tions for the macroscopic amplitudes up to an arbitrarily high order of the scaling 
parameter, thereby deducing the resonance and non-resonance conditions on the 
fixed wave vectors and frequencies of the pulses, which are required for that. The 
derived equations are justified rigorously in time intervals of macroscopic length. 
Finally, for sets of up to three pulses we present a complete list of all possible 
interactions and discuss their ramifications for the corresponding, explicitly given 
macroscopic systems. 

Key words and phrases: nonlinear discrete lattices; interaction of modulated pulses; 
multiscale ansatz; derivation and justification of macroscopic dynamics. 
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1 Introduction 

In recent years a big part of the research activity within applied mathematics has been 
focused at the investigation of so called multiscale problems, cf. for a survey [MieOGb] . One 
aspect of this area is to investigate for a given, mostly physical, phenomenon the relation 
between its descriptions at different time and/or space scales. A prominent paradigm 
concerns the dynamics in lattices. Here, the underlying microscopic model is a lattice of 
atoms interacting with each other, subject to some given potential, and possibly embedded 
in a background field. Then, one has a complete description of the dynamics within the 
lattice by the discrete system of Newton's equations of motion for each atom (see (11. 2p ). 
However, thinking of the distance between the atoms as very small, one is interested in 
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the dynamics of some macroscopic object within the lattice. For instance one could ask 
how initial data varying on a large space scale as compared to the distance of the atoms 
evolve in the lattice. The corresponding dynamics would be described by a macroscopic 
model, which in the limit of the distance between the atoms going to zero is a continuum 
system. 

Of course, this idea is known within the physics community since long, and macro- 
scopic, or effective, continuum dynamics have been formally derived from microscopic, 
discrete lattice systems for a plethora of phenomena, usually by some sort of perturbation, 
or asymptotic, expansion with respect to the scaling parameter < e -C 1 of the distance 
between the atoms. 

However, the derived models mostly lacked a justification in a mathematical-analytical 
rigorous sense. Moreover, at latest since the unexpected numerical discoveries of Fermi, 
Pasta and Ulam in 1955 (see |FPU55j ) on the behaviour of waves in nonlinear oscillator 
or atomic chains, i.e., one-dimensional lattices, it was clear that the understanding of the 
dynamical behaviour of even one-dimensional lattices is far from complete. This surely 
played a role in arising the interest also of the mathematical community and leaded to an 
increased research activity regarding the dynamics in discrete lattices. We refrain here 
from giving an overview over the huge amount of work done since then, and in particular 
the various lattices and potentials considered, the (macroscopic) dynamics derived, and 
the mathematical techniques used, for which we refer the reader to the surveys |CRZ05l 
IDHROGt IGHMOGj and the references given therein. 

Instead, we present in the following precisely the lattice model considered, the ques- 
tion discussed, the method used, and the result obtained in the present paper, and refer 
later on to literature related to them. 

The nonlinear lattice. We consider the d-dimensional (Bravais) lattice F C M'^, d G N, 
r = {7 = ^1^71 + ■ ■ ■ + kdOd : k = [h, ...,kd)e U"}, (1.1) 

where {gi, . . . ,gd} is a set of linearly independent vectors of (yielding F = Z'^ when 
{9i)j = ^ij)^ the microscopic model 

= J2K{^i+o.{t)-x^{t)) - W'{x^{t)), 7 e F, t G M. (1.2) 

This is a system of infinitely (countably) many coupled second-order ordinary differential 
equations, describing according to Newton's law the scalar displacement a:^(t) G M at time 
if: G M of a particle (e.g. an atom) of unit mass 1 from its rest position 7 G F due to the 
(pair-)interaction and on-site potentials : M ^ M, a G F, and : M ^ M, respectively. 
One can think of each mass particle of the lattice of rest position 7 as being connected to 
every other particle 7+a, a G F \ {0}, via a spring of force V^, and of the whole lattice 
as being embedded in a background (external) force field W. Since the interaction and 
on-site potentials Va and W are independent of the position of the particles, and since 
all particles have equal mass, we speak of a monoatomic, homogeneous lattice. Moreover, 
since naturally V^{x^+a—x-y) = —Vi^{xy—Xy+a), we assume V^(x) = V-a{—x) for all 
a G F, X G M. Note that we allow for interactions among pairs of particles, which are 
arbitrarily far apart from each other. In order for our results to hold true, we will need 
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to impose some rather strong decay conditions (cf. (13.71) ) on the interaction potentials 



VJx 



as \a\ 



oo, which however seem physically plausible. Of course, these 
conditions are satisfied if we assume that pairs of particles interact only for distances up 
to some finite range i? > 0, i.e., Vq, = for |a| > i?. For instance, for F = Z and = 
if |a| > 1 we obtain the classical oscillator chain with only nearest-neighbour interactions 

x^{t) = Vl{x^+i{t)-x^{t)) - Vl{x^{t)~x^_a{t)) - W'{x^{t)), 7 G Z, teR, (1.3) 

which for W = and anharmonic Vi, Vi(x) 7^ 0-1,1^^ is the Fermi-Pasta-Ulam (FPU) 
chain, and for harmonic Vi and anharmonic W the so called Klein-Gordon (KG) chain. 

Throughout this paper we assume Va, W G C^''"^(M) for some G N to be specified, 
and V^(0) = V^(0) = W{Q) = W'{0) = 0, which allows for the Taylor expansions 



N 



V;,{x) = Y,(^n,aX'' + V^ 
n=l 

N 

W'{x) = ^6„x" + iy; 
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^Af+l,o(^) 
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iN+2) 



X] 



{N+l)\ 

W^^+^^{kx) 



-X 



N+l 
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-X 



N+l 



;i-4) 



with an,a := V^i""*''H0)/(n!), 6„ := l^("+i)(0)/(n!) and ^ = fe(x), k = K^ix) G (0,1). 
Note, that Va{x) = V-a{—x) implies a„,Q = (— l)"+^a„_a. By (11.41) . the microscopic 
model (11. 2p reads equivalently 
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5^ 5^ an, 

n=i Voer 



^nX^ 



W^^,{x^] 
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Finally, note that (11.21) is an infinite-dimensional Hamiltonian system with conserved 
total energy (Hamiltonian) 7i(x, x) = /C(x) + U{x), and kinetic and potential energies 



)C{x) 



7er 



U(x) 



7er \ aer / 
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Pulses and their interactions. The linearized microscopic model (11.21) 



possesses the plane-wave solutions or pulses x^(t) = E(t,7) + c.c, where c.c. denotes the 
complex conjugate of E(t, 7) := Q^i'^^+'^-y) ^ provided the frequency a; G M and the wave 
vector -i? G Tr satisfy the dispersion relation 



00^ 



:= (l-e^'-") +bi = Yl (1- cos(?9-a)) + h (1.7) 



(the latter equality follows from above). Here, Tr := M!^/r* is the d- 

dimensional (dual) torus associated to the lattice F, where = Lin(]R'^) is the dual of 
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M'^' and r, := {^9 G : d : a ^ 'd ■ a, 'd ■ a e 2tiZ V « G T} is the dual lattice to T. In 
the following we assume that the stability condition 



> for all 7? G Tr 



(1.8) 



is satisfied, and set the dispersion function > for all "(9 G Tp. Note, that (11.81) 

implies that there are no pulses with frequency u = 0, and in particular that E(t,7) = 1 
is not a pulse. The stability condition can be guaranteed by appropriate choice of the 
coefficients ai^a and bi of the harmonic parts of the potentials Va and W. For instance, 
it is surely satisfied if we assume min{6i, 2 y minjai g, 0} + bi} > 0. 



The dispersion relation enables the following characterization of pulses: A pair {-d, u) G 
Tr X M represents a pulse if and only if 6 := — cj^ = 0. The pair — ("i?, uj) represents 

the same pulse, while ±(-(9, —u) represent the associated pulse traveling in opposite direc- 
tion. In the following we will denote pulses mostly by the pairs representing them. In fact, 
generalizing the use of this notation, we will denote arbitrary functions E(t,7) = e'^'^*"*"''''^-' 
by the pair ("i?, o;) G Tr x R. Thus, be aware that this notation does not imply that E is 
a pulse, unless the corresponding pair satisfies the dispersion relation (11.71) . 

In this paper we will deal with sets of z/ G N different, fixed pulses {(-i^j, cu-,) : j = 
1, . . . , I/}, such that 5j := n'^i^j) - uj] = and E^ := e^K'+'^^T) ^ E^ for j ^ i. With 
the notation {'d_j,uj^j) := —{{}j,Uj) this set of pulses can be written equivalently as 
{{-djjUj) : j G A^}, where Af := {— z^, . . . , —1, 1, • • • , i^}, and it holds E_j = Ej, z denoting 
the complex conjugate (c.c.) of 2; G C. 

Since the microscopic model (II. 2p is nonlinear, products of pulses E(jj^...jj.) := Ej^- ■ -Ej^ 
with (ji, . . . ,jk) G JV^ will play an essential role throughout this paper. However, different 
index-vectors (ji, . . . can yield the same product. (Confer for instance the products 
corresponding to the indices (p, g), {q,p), or {p,—p,q), {—p,p,q), {p,(l,—p), {—p,(l,p), 
{q,p,—p), {q,—p,p) for p,q E Af.) Moreover, the product for (ji,...,jfc) can equal a 
product for (ii, ...,««;) with n ^ k. (Confer for instance the products corresponding to 
the indices (p, — p, q) and q.) Hence, in order to identify each appearing product of pulses 
by a unique index, we introduce the following 

Notation 1.1 

(a) We denote by Jm a representant of all those (ji, . . . , j^) G Af'^, k < k, which lead 
to the same product Ej^{t,'j) := e'^'^-^^^^^^^rn-i) ■= E(j^^...j^)(t, 7), where § : = 
+ . . . in Tr and uj^ := ujjj^ + . . . +ujj^. The set of all (different) represent- 
ants Jm within |Jk=i is denoted by T^. The index m denotes Jm £ Af'^\Tm-i with 



(b) The representation Ej-^ = E(j^ . j^) is denoted by Jm = [ji, ■ ■ ■ , Jk)- More generally, 
= ^ih,...,i.) denoted by (ji, . . . , j^) = (i^, . . . , v). 

(According to this notation, the products E(p = E(gp) from above have the same repre- 
sentant Jm = {p, q) = ((?, p) with m = 1 iff (p, g) = j G and m = 2 iff (p, q) 7^ j for all 
i G AT. Moreover, E(p _p_q) = Eg has the representant J7i = g G A/", as do all permutations 
of the indices within (p, — p, q).) 

By use of Notation II. H we can state concisely some observations on sets of pulses. 




Ji:=jGA' = Ti. 
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Remark 1.2 Let {(-i^j, cUj) : j G N"}, M = {±1, • • • , ±z^}, be a set of v pulses. Then 

1. It always holds U^'"^ C , since (ji, . . . , jfc-2) = (ji, • • -JkJ, -j) for any j e U . 
In particular, (ji, . . . , jj^) G can have a representant Jm. G Tk-i- (Cf., e.g., for 
k = 3 the example {p, —p, q) above.) 

2. In general M^-^ (JL . E.g., for N = {±1}, we get Af^ = {±(1, 1), ±(1, -1)} and 
%nM'^ = {±1, ±(1, 1, 1)}. Since E±i ^ ^ E^^ and E±i EiE_i = 1 ^ E^^, 
we obtain A/" n A/"^ = and Af^ nAf^ = 0, while Af C Af^, in line with Obs. [H 

3. If for each j E Af there exist p,q E Af such that {p, q) = j, then Af^^^ C Af'', since 
(ji, • • • , jfc-2, j) = (ji, • • ■,Jk~2,P, q)- In this case Tk C Af''. 

In general, a product of k pulses E^-^ . . . E^,, need not be a pulse. However, if it is a 
pulse, we make the following definition. 

Definition 1.3 A; G N pulses {-t^j-jLOj-) , i = l,...,k, are in resonance (of order k) or 
interacting with each other (in order k) if 

or, equivalently, by use of Notation \l.l\ if 

5^™:=^^'(^^,J-4„=0, (1.9) 

where Jm = (ji, ■■■,jk)- 

Finally, we conclude our discussion of pulses, their products and their interactions, by 
providing a definition which will be essential for our results. 

Definition 1.4 The set of pulses {{i^jjUjj) : j G Af} is closed under interactions up to 
order k (or for short: fc-closed ) if 

h^--=^'i^jJ-^L^O forJme%\Af. (1.10) 

In other words, a set of pulses satisfies the closedness condition f ll.lOp if none of the prod- 
ucts of up to k (possibly identical) pulses of the set are new pulses, i.e., pulses which are 
not included in the set. As the negation of (11.91) . the conditions (ll.lOp are refered to also 
as nonresonance conditions. 

Macroscopic dynamics of modulated pulses - Formal derivation. For a given 
pulse E(t, 7) = e^^'^*"'"'' '''^ an (amplitude-) modulated pulse is a function of the form 

x^{t) = e'' A{eH, e{-f-ct))E{t, -f) + c.c. 

with the scaling parameter < e <^ 1 and the amplitude or envelope {e°')A : R x M'^ — » C. 
The latter depends on the macroscopic space variable y = 5(7 — at), where we allow for 
possibly moving-frame space coordinates of velocity c G M'^. The dependance of A on 
such a y means that the amplitude varies at a much larger space scale of order 0{-) — 
or, for space-periodic A, that it has a much greater wavelength — in comparison to the 
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space (wavelength) scale of order of E(0, 7) = e^'^"^, since < e <^ 1. The choice of 
the powers a and b determines the size of the amplitude and the length of its time scale 
r = e^t. This choice is motivated by the physical phenomena one wants to study, and 
shapes the macroscopic equations obtained, cf. the discussion of related literature below. 

Here, for a set of G N different, fixed pulses : j e A/"} we consider their 

small-amplitude macroscopic modulations 

V 

iX^%{t) :=£ J]A,(£t,£7)E,(t,7) + c.c. (1.11) 

with hyperbolically scaled amplitudes Aj : RxM*^ — > C: t = et, y = ej. We are interested 
in the macroscopic dynamics of such modulated pulses in (11.21) . This means that, making 
for solutions to the microscopic model (II. 2p the multiscale ansatz 

x,it) = {xt%it) + Oie'), (1.12) 

we are interested in the dynamics of the amplitudes Aj, or more generally, making the 
ansatz 

N 

x,{t) = {Xi%{t) + 0{e''+'), Xi'^:=j2e' ^ A,,j^Ej^ (1.13) 

k = l JmGTfc 

for G N, in the macroscopic dynamics of the functions Af^^j^ : MxM'' C, (r, y) 
Ak,jmi^,y), (with Ak,j^ = Ak-j^ in order to obtain X^'^ G M), and Aij = Aj for j eM 
(recall the definition of in Notation ll.il) . Note here, that in view of Remark 11.21 Obs. 
[21 for each k we take into account Jm G instead of Jm £ (cf. also Sec. 14.1.31) . 

Inserting (11. 13^ into the microscopic model (11. 2p and carrying out the operations, we 
obtain on the left- and right-hand sides expansions of the form (11.131) . where, in particu- 
lar due to the nonlinearity of (11.21) . the macroscopic coefficients are products of functions 

Jm , < their space and time derivatives. Thus, in order for the microscopic system 
(II. 2p to be satisfied up to order by the approximation (11.131) . and since at each k the 
microscopic patterns E are different from each other for G and are nonvanishing 
(of modulus 1), the macroscopic coefficients on the left- and right-hand sides of the cor- 
responding microscopic patterns have necessarily to be equal. This yiels an hierarchy of 
equations for the functions Ak^j^ in (I1.13p . which guarantee that X^'^ satisfies (II. 2p up to 
order e^, and constitutes the formal derivation of their macroscopic dynamics, carried out 
in Section [2l In particular, due to the hyperbolic scaling t = et, y = e'y and the scaling 
by 6 of the first-order amplitudes Aj, it turns out that the latter are determined by the 
equations for e'^Ej, which contain the macroscopic first-order time- and space-derivatives 
Ot-Aj and VyAj, and, due to the cubic terms in the potentials Va, W, the products ApAg 
for which EpEg = E^, whenever the pulses {{}p,Up) and {{}g,ujg) interact to generate the 
pulse {{}p,Up) + {'dq^ujq) = cf. (I2.17p . Note here, that due to the underlying hy- 

perbolic scaling and the scaling of the amplitudes by e, coupling of more than two pulses 
is not shown in the leading order effective dynamics. Moreover, as we will see in more 
detail in Section [21 in order to calculate the functions A^^j^, k < N, G T/v \ A/", we 
have to require that the system is closed under interactions up to order A^. 
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Justification. The macroscopic evolution equations we obtained by tlie above procedure 
of tlie formal derivation, however, establish only the necessary conditions on A^^j^, which 
solutions to the microscopic model f 1 1.2 1) have to satisfy z/they are of the form fll.l3p . The 
existence of such solutions is not at all guaranteed. Indeed, there exist counterexamples, 
cf. jSdl95l[Sdl05] . 

Since for t = we can prescribe the initial data in such a way that they have the 
form fll.l3p . the question of the justification (or the validity) of the formally derived 
macroscopic equations is answered in the affirmative, if we can show that solutions to the 
microscopic model (11. 2p with such initial data maintain this form also on time intervals 
[0, To] of positive macroscopic length tq > 0. The reason why such a condition should 
be satisfied on the macroscopic time scale r = et is obvious, considering that the scaling 
parameter is very small, < e <^ 1. 

For the systems considered in the present paper, the justification of the macroscopic 
equations derived formally in Section [2], is carried out in Section [31 In particular, the 
justification result. Theorem 13. 6[ states that if the error between the approximation X^'^^ 
in (I1.13P obtained by the solutions A^^j^ to the macroscopic equations (viz. (I2.17p . (12.180 
for N = 2 and (I2.20p . (I2.2ip for > 3) and a solution of the microscopic model (11.21) is 
initially of order with (3 G {1, N — d/2], then it will remain in the same order for times 
t < To/e, Tq > restricted only by the existence of solutions of the macroscopic equations. 
In particular this is the order of the error when the initial data for the microscopic system 
(ll.2p are given by ^^'f i- Thus, naturally, the higher the order of the approximation the 
smaller its error with respect to a true solution. However, note that there is a lower bound 
for the order N in the ansatz (I1.13P needed in order to get reasonable results, depending 
on the space dimension d of the lattice and the scaling by e of the amplitude. Hence, 
only in the case of a one-dimensional lattice, i.e. an oscillator chain, the order N = 2 
is sufficient to obtain valid macroscopic dynamics of the first-order amplitudes. For the 
other physically relevant space dimensions = 2 or = 3 we need to use = 3 in (11.130 
in order to be able to determine the full second-order approximation X^'^ which gives the 
valid macroscopic dynamics. 

The proof of Theorem 13.61 consists on a Gronwall argument for the error = e~^{x — 
X^'^) between approximation and original solution x of (II. 2p in the energy norm of the 
microscopic system. This needs, on the one hand, an estimate of the residual terms 
res(X^'^) = 0{e'^~^^) in the in the £^(r)-norm. Since they depend on the solutions A^^j^ 
of the macroscopic equations, we obtain due to rescaling ||res(X^'^) ||£2 = 0{e^~^^~'^^'^) 
provided the regularity of A^^j^ is sufficiently high. On the other hand, we have to 
estimate the difference between the nonlinear terms of the original solution x and of the 
approximation X^'^ in the integral formulation of the differential equation for i?^, see 
(I3.19p . Since they are quadratic, they can be matched by the hyperbolic scaling of the 
system where r = et, which is essential for the Gronwall argument, the latter yielding 
finally the justification of the macroscopic equations for r < tq. 

The strategy of the proof is classical within the theory of modulation equations and 
not restricted to discrete systems. For a concise presentation in the case of continuum 
systems, see |KSM92] . For a further example on how essential the matching between 
macroscopic time scale and nonlinearity is, cf . |GM04l IGM06] , where the dispersive scal- 
ing r = eH, y = e{-f - n'{^)t) was used in order to derive the nonlinear Schrodinger 
equation for an oscillator chain. As long as the nonlinearities of the system are cubic. 
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[GM04j . the proof is straightforward. In order to include cubic potentials, [GM06j . one 
has first to apply a normal-form transformation on the system. In the present case, due 
to the hyperbolic scaling, the latter is not needed. 

Examples. We conclude our paper with a list of examples in Section |H In particular, 
we classify all possible interactions (or order 2) for up to z/ = 3 pulses which form systems 
closed under interactions up to order N = 2, and give the corresponding macroscopic 
dynamics. The main purpose of the section is on the one hand to show how the resonances 
are mirrored in the macroscopic dynamics and on the other hand to clarify the role of the 
closedness condition. The main observation is that if a system is not closed, i.e. if the 
considered pulses generate new pulses not taken into account by the multiple scale ansatz 
f ll.lSp then the macroscopic equations describe only the dynamics in trivial cases. In other 
words, in order to be able to detect the dynamics of interacting pulses macroscopically 
one has to take into account all pulses generated by the considered system. However, 
generating pulses can of course be ignored. 

As explained above, for one-dimensional lattices the requirement on the systems of 
pulses to be closed under interactions up to order = 2 is sufficient in order to obtain 
valid macroscopic dynamics. For two- and three-dimensional lattices the systems need to 
be closed under interactions up to order = 3. This increases the number of nonres- 
onance conditions fll.lOp which have to be satisfied in order to obtain the second-order 
approximation X^'"^ needed for the justification result. However, it leaves the resonance 
conditions and thus the macroscopic equations for the first-order amplitudes unchanged. 
Exemplarily, we give the macroscopic equations for the second-order amplitudes A2J 
J = 1, . . . , 1/ for a single pulse (Sec. I4.1.3P and for the three-wave-interaction (Sec. I4.3.3p . 

Finally, in Section 14.41 we look closer at the resonance condition for a three-wave- 
interaction in a one-dimensional lattice with only nearest-neighbour interaction and a 
stabilizing on-site potential. In particular, we determine the coefficients of the harmonic 
part of the interaction potential which allow for the existence of three-wave-interactions. 
It turns out that indeed this is possible only for a small range of repulsive harmonic parts. 

For a more precise discussion and interpretation of the phenomena observed in the 
interaction of pulses, see the introduction to Sec. HI 

Related literature. The present work is closest related to work which also uses the 
modulational approach described above in order to justify modulation equations derived 
formally by asymptotic expansions in the scaling parameter. For a general description of 
the method in the case of lattices, see |GHM06] . As already mentioned, this method was 
used also in |GM04tlGM06j . justifying the nonlinear Schrodinger (nlS) equation describing 
for a single pulse in an oscillator chain the deformation of its (small) amplitude in the 
macroscopic variables r = e'^t, y = e{'j — fl'{i))t), the microscopic group velocity 

of the carrier plane wave. In this case the dispersive scaling, i.e. in particular the longer 
time scale was chosen in order to allow for the amplitude scaled be e (and 

corresponding to a weak nonlinearity) to deform, and as a consequence the nlS equation 
was justified for t G [0,ro/e^]. Under the hyperbolic scaling used here, we observe only 
the transport of the amplitude, see fl4.2p . 

As indicated previously, the modulational approach is not at all restricted to lattice 
systems. In the contrary, it has been previously applied to continuous systems, see for an 
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overview |Kal89t ISchOS] and for a short exposition of the idea |KSM92j . Concerning the 
interaction of pulses, in [SW03] it is shown that resonating water waves which are subject 
to weak surface tension can be approximated by using the same scahng as in fll.l3p . by 
a system of three- wave-interaction equations of the form (14. 141) . For a discussion of this 
system and its various physical applications, see |BS90t IKauSOl IKea99l IKRB791 ISW03j 
and the references given therein. Moreover, we would like to point out the structural 
similarities of the derived three- wave interaction equation fl4.14l) and its underlying setting 
with the Boltzmann-like equation describing the collision of phonons in the kinetic limit 
and its relations to wave turbulence, see |Spo06| . 

In |SW00j a coupled system of Korteweg-de Vries equations was justified as governing 
the evolution of the amplitudes / and g of two counter-propagating waves in an FPU 
chain in the "small, long- wavelength limit" (cf. also |FP99] ). i.e., making the ansatz 

x^{t) = eV(^(7 + ct),eH) + e^g{e{-f - ct),eh) + 0{e^). 

While the previous ansatz contains no internal microsctructure, in |SUW09j the interac- 
tion of two (weakly) amphtude modulated pulses of different group velocities Ci 7^ C2 and 
time-independent amplitudes is considered in the chain fll.3p under the dispersive scalings 
T = e^t y = e{'~f — Ci^2t), and it is proved that after interaction the amplitudes retain their 
shape but experience a shift of order 0{e) in position; see also |BF06] . Analogous results 
are obtained in |CBea07[ ICBeaOSj on a continuous one-dimensional string. 

In a different physical setting but with the closely related WKB-approximation ap- 
proach, in [GMSOSj the amplitude equations for interacting modulated pulses of a nonlin- 
ear Schrodinger equation with periodic potential are justified in the semiclassical scaling. 

The results obtained in the present paper rely for the derivation of the macroscopic 
equations fully on Newton's equations of motion (11.21) . However, this microscopic model 
possesses Lagrangian and Hamiltonian structure (LHS), cf. (11.61) . In |GHM08a] a general 
framework is presented for the direct reduction from the microscopic LHS of the macro- 
scopic LHS corresponding to the limit equation by use of the related two-scale ansatz. 
Among several examples (with or without microstructure) this is also performed for the 
three- wave interaction of Sec. 14.3.11 Case [2] (see also [GHMOSb ]). 

We conclude our (non-exhaustive) survey of related literature by recalling three fea- 
tures of the setting we discuss in the following: first, the assumption of scalar displacement 

G M, second, the existence of a stabilizing on-site potential, and, third, the smallness of 
the amplitudes (weakly-nonlinear regime). These features imply also possible generaliza- 
tions. In this direction, |Mie06a] contains a thorough analysis of macroscopic continuum 



limits in multidimensional linear lattices, while [DHM061 IDHR06] present first results con 



cerning the Whitham modulation equation as the macroscopic limit in oscillator chains 
in the fully nonlinear case. 



The paper is organized as indicated by the introduction: in the following Section [2] 
we derive formally the macroscopic equations for the functions Ak^j^ of a general ap- 
proximation (11.131) . which are then justified in Section [31 In Section H] we present all 
possible macroscopic systems (mainly) for the first-order approximations for systems of 
up to z/ = 3 pulses, discuss exemplarily the significance of the closedness condition, and 
prove in a typical case the existence of interacting pulses in oscillator chains. The main 
observations are summarized in the beginning of the section. 
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2 Formal derivation 



In this section we derive the macroscopic equations which the functions Afc,j,„ of the 
multiscale ansatz X^'^, see fll.lSp . have to satisfy necessarily in order for the microscopic 
model (11.21) or, equivalently, (11.51) to be satisfied up to terms of order e'^ . We follow the 
procedure outlined in the Introduction. We chose to present the derivation in full detail, 
in order to enable the interested reader to see clearly the origin and structure of the 
resulting macroscopic equations in their general form. However, the reader interested only 
in the equations themselves may proceed directly to equation (12.121) . which summarizes 
the obtained results. We present the macroscopic equations explicitly in order to provide 
them for any possible further use, as for instance for the derivation of the macroscopic 
equations in the concrete examples of Section HI 

We start by inserting the ansatz X^'^, see (11.131) . into the microscopic model (11.51) . 

N 

n=l aer 

where (x.+q)^ := x^+q for x G ^^(r). Next, we expand the left- and right-hand sides of 
(12. ip with respect to terms of the form e^Ej^, k = 1, . . . , N, e (cf. Notation [LTD- 
For the left-hand side one has 

N 

N+2 

= ^£%, h:= {d'rAk-2,j^ + 2iujj^drA,_^,j^-uj'j^A,,jjEj^, (2.2) 

k = l Jn,&Tk 

where A^^j^ = for k G Z\{1, . . . , A^} or m > k. For the expansion of the right-hand side 
of we first expand (X^'^).+„-X^•^ Since {Ak,j^EjJ.+^=Ak,jJ; ■+ea)e'^^^-^Ej^, 
we use the Taylor expansion of Afc,j™, with respect to y, assuming for the moment that 
Ak,jjr, ■) G C^-^'+HM'^) for r G [o',ro], tq > 0. Hence, we obtain 

N 

ixin-+.-Xi'^ = J2e' Yl iAk,j,A;-+ea)e^'--"-A,,j^)Ej^ 

k = l Jm&Tk 

N /N-k+1 \ 

fc=i JmeTk \ S=0 J 

N+1 k 

= E^^'^'''"' dk,a ■= 'Y E '^^'*=~^-^'"'"Ejm; (2-3) 

k=l Jm&Tk l=m 
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with 



D 



s\ 



s = 0, 



s = 1, . . . , A^-fc, A; < A^, 



(2.4) 



^AT-fc + l^Cl) 



s = N-k+1, 



where := Ak,j^{-,-+hea), h = hk^j^^e^{ea) G (0,1). (Note, that L'jv+i,o,j„,a = 

and Dk s,j^,a = for m > k.) Next, abbreviating 



N 



X 



A,e 
N 



(2.5) 



(cf. (11.131) ). we use (12.31) . (12.51) in order to expand the first sum on the right-hand side of 
(12. ip in terms of ^ 



N 



'iV+1 



N 



n=i Vaer \fc=l / \fc=l 

AT n(Af+l) / n n 



au. 



n=l k=n fci + ...+fe„=fe \Q,gr j=l 

l<fc,-<JV+l 



i=l 



N{N+1) 

E 

fc=i 



k 



min{k,N} / n n 

Y Yi XI n ^-^-^ ~ n ^'^^ 



[ fc-i 1 I -I fei+...+fcn=fc \aer 2=1 
Ljv+iJ^ i<fcj<]v+i 



2=1 



(with X = [x] + /i, [x] e Z, /i G [0, 1) for x G M). It remains to obtain the expansion of 
in terms of Ej-^. By (12.51) . (12.31) we get for the summand in brackets in Sk 



Jrr 



agr i=l Jm^eTk^ £i—mi 



i=l Jm.GTfc. 



\ 



ei+...+en<k a&r i=l 



i=l 



and hence by inserting the last expression in the definition of Sk above 

min{fc,Af} 

E E E E 

n n 
1=1 ei+...+i„<k-i agr i=l 
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with 

n 

cu.„...,^.J :=E«"."n(^''"'"^'"-l) (2-7) 

aer i=i 

In particular for n = 2, recalling 02,0 = — a2,-« for a G F, we obtain 

C{j^,J^.) = -4i 2^ a2,a sm (^^-a J sm J sm ^ a j - 62- (2.8) 



For k = 1, . . . ,N the involved expansion fl2.6p can be structured by separating the terms 
involving the functions Ak^j^ and Ak-i^j^ from all the others, which we subsum into the 
term Sk~2,jm- This structure is essential for the derivation of the macroscopic equations 
for j^. Indeed, by writing out the summands for n = 1 and n = 2 in (I2.6P when 



k = 1, . . . , N, we obtain 



+ min{A;-l,2} ^ c^p,j,)Ai,pAk-,,j^ + Sk-2,j^'^Ej^ (2.9) 



with 



fc-2 k-2 



Sk~2,J^ X] X] + ^ (^{J^^,J^^)AK,J,,,^Ak-n,J„, 

i=m K=2 (Jmi.Jm2) = Jm 

2 

+ Y Xl"2,an^^-'=^-^-^™>'" 



fcl+fe2 = fe (Jmi,Jm2)^J'm ll+e2<k~l agF 1 = 1 

k 



l<ki<k-l Jm.i'^T'ki ^i<h<l^i 



EE E 

n=3 fcl+-.. + fcn = fc {Jm-^ Jmn) = Jm 

l<fci<fc-n + l Jm.'iTu. 

C{Jm^,-,Jm„) W Ak^^J^^ ^ Y Y n ^f-^,k,-t,,Jm.^,a ■ (2.10) 



i=l ii+...+ln<k-i a&T i=l 



For ([23]) we used 

$^ai,„Z}fc_i,i,^„,, = Y,ai,J'>^--^a-VyAk.,,j^ = 2in{d jjV^^j J -VyAk-.^j^ 
aer aer 

(2.11) 

(cf. (12. 4p and (11.70 ) for the term obtained for n = 1, and the formula 

k-l k-2 

Y^ ^f^^k-K = min{fc-l, 2}ai,fc_i + Y (^>i,k-K for ^K,k — K ^k — K,K 

K=l K=2 
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for the term obtained for n = 2. (Note, that Sk-2,jm contains indeed only functions A^^j^ 
with i<k-2.) Finally, since V^^^^Jx) = 0{\x\^+^), W^^^{x) = 0{\x\^+^) (cf. nl 
and (X^'").+„-X^'" = 0{e), X^'' = 0{e) pointwise for (^,7) e [0, cx)) x T (cf. (IXSD . 
f ll.lSp ). the second sum on the right-hand side of (12.11) is of order 0{e^~^^). 

Using the above expansions, (12.11) reads 

AT 

Y^e'^itk - Sk) + res(X^'^) = (2.12) 

k=l 

with tk, Sk given by (12. 2p . (12. 9p . and the residuum res(X^''^) = 0{e^~^^) given by 

N{N+1) 

res(XiJ'^):=e^+W + ^^+W- ^ '"'^ 

k=N+l 

-Y.^Ua{iXin-^c.-Xi'^) + W'^_,,{X^'^) (2.13) 

with Sk given by (12.61) . Thus, the approximation X^'^ satisfies the lattice system (11.51) 
pointwise for {t, 7) G [0, 00) x T up to order if and only if 

4 - Sfc = for all k = l,..., N. (2.14) 

Furthermore, by (12.21) and (12.91) . tk — Sk is given for each A; = 1, . . . , iV as an expansion in 
terms of the harmonic functions Ej^, J7m G Tk. Since the coefficients of this expansion, 
the macroscopic amplitudes Ak^j^ and their derivatives, are varying much slower in space 
and time as compared to the microscopically oscillating, mutually different, non- vanishing 
patterns Ej-^ (with different wave- vectors and frequencies and modulus 1), in order for 
the expansion tk — Sk to equal 0, each of its coefficients has to vanish identically. Hence, 
the approximation X^'^ satisfies the lattice system (11.50 up to order if and only if 

hrr.Ak,j^ + 2iuj^d^Ak^,,j^ - 2in{^jjVM^jJ-'^yAk-i,j^ 

- mm{k-l, 2} ^ c^p^j^)Ai^pAk-i,j^ = Sk^2,j^. - d^Ak-2,jr^ 

for all = 1, . . . , X and all Jm G Tk, (2.15) 

where we used -u}^^^ - c^j^) = -u}^ + n'{i!}jj = 6j^, cf. (EZD, dEZD and The 
equations (I2.15P are the macroscopic equations, which the functions Ak^j^ of the approx- 
imation X^'^ given by (I1.13P have to satisfy necessarily, in order for this approximation 
to solve the microscopic model (11.51) up to residual terms of order e^"*"^. However, the 
equation for Ak-j^ is just the complex conjugate of the equation for Ak^j^. Thus, since 
Ak-j^ = Ak^j,^, it suffices to determine just one of them. 

Since the equations for any given k = 1, . . . , N involve the functions A^^j^ with k < k 
it is natural that they have to be solved inductively for increasing k. Hence, recalling 
that Ak^j^ = for G Z \ {1, . . . , N} (and thus in particular Sk-2,jm = for /c < 2) and 
m > k, the equations (12.151) for k = 1 read 

5jAij=0 for J G A/". (2.16) 
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Since by assumption Sj = — C(j| = 0, reflecting the fact that the flrst-order approx- 

imation X/^'^ in (11.131) consists of modulated pulses, these equations are automatically 
fulfllled and the amplitudes Aij remain at this stage undetermined. Of course, these 
equations can be interpreted also in the opposite direction, namely as requiring necessar- 
ily from the pairs ±{'dj,u!j) to satisfy the dispersion relations Sj = 0, i.e. to characterize 
pulses, in order to allow for non-vanishing amplitudes Aij in the approximation . 
Then, the functions Aij are determined by the equations fl2.15p for A; = 2 with 6j = 

2iujjdrAij -2in{^j)V^n{^j)-VyAij = ^ C(p,g)Ai,pAi,g for j G AT (2.17) 

p,geJ\f 



and 



'^^2^2,^2= hi J2e^^''\^^, (2.18) 



where by (12. 8p 




(2.19) 



The equations (I2.17P determine the amplitudes Aij, j G Af. If for some j E Af there 
exist p,q E Af such that (p, q) = j then the corresponding equations are semilinearly 
coupled transport equations. If for some j E Af there do not exist such p,q E Af, then 
the corresponding equation for Aij is just a transport equation with vanishing right hand 
side, uncoupled from the other equations. The same applies if the coupling coefficient C(p g) 
given by (I2.19p vanishes, for instance when the interaction and on-site potentials Vq, and 
W, cf. (II. 4p . have no cubic terms. Of course, within the considered system of pulses Af 
there can exist subsystems of pulses interacting with each other but not interacting with 
other (subsystems of) pulses. The corresponding macroscopic equations then establish 
(closed) coupled subsystems, cf. also Sec. HI Note that each amplitude Aij is transported 
by the group velocity V^fl{i!}j) of its corresponding pulse. 

Having determined the amplitudes Aij for j G A/" by (12.170 . we are then able to 
calculate the functions ^2,^72 <^2 G A/'^\A/'by (12.180 . provided 7^ 0. If Sj^ = 0, then 
ii^'d j^) = ±{^p + i!^q,ujp + ujg) characterizes a pulse not considered (or, equivalently, 
assumed to have an identically vanishing amplitude) in our approximation Xf'"^, and 
(I2.18P can be seen as a further condition on the first order amplitudes Ai,p, Ai q for which 
{p, q) = J2- As we will exemplify in Section HI this condition then implies that some of the 
involved first order amplitudes have also to vanish identically, thus effectively prohibiting 
the description of the macroscopic dynamics of the corresponding pulses. This problem 
can be overcome if we include into N all (p, g) G A/"^ with 5(p y) = 0. Then ^ Q 
for J2 G A/"^ \Af, and we call such a set of pulses {±{'dj,ujj) : j = l,...,z/} closed 
under interactions up to order k = 2, according to Definition 11.41 Then, for a set Af of 
pulses with this property, we can calculate by (I2.15P for = 2 the functions ^2,^2 
J^2 € Af^ \f^y while the second order amplitudes A2J, j G Af, remain undetermined by 
these equations. 



14 



The appearing pattern, namely that the functions A^^j^ with J'^ E Tk\Af can be 
determined by the equations fl2.15l) for k and the k-th order amphtudes A^j with j E M 
by the equations for /c + 1, can be continued inductively, provided the set of pulses N 
is closed under interactions up to order k. More precisely, for k = 3, . . . ,N, and using 
6j = for j e J\f, the equations (I2.15p read 

2iujjdrAk^ij - 2in{'dj)V^n{'dj)-VyAk-ij - 2 ^ C(^p^q)Ai^pAk-i,g 

= 2 c^p,j,)Ai,pAk-i,j^ + Sk-2,j - dlAk-2,j for j eN (2.20) 



and 



h^A.^Jm = -'^^^j^drAk-i,j^ + 2in{^jjVM^jJ-'^yAk-i,j^ 

+ 2 Yl c^p,j,)A,,pAk.i,j^ + Sk_2,j^ - d^^Ak-2,j^ ioiJme%\M. (2.21) 

The equations fl2.20p form a system of (in general) linearly coupled inhomogeneous linear 
transport equations for the (k — l)-th order amplitudes A^-ij, j G A/", travelling again 
with the group velocity of the pulse they modulate. As in the case for k = 2, the equations 
for j E M for which no p, g G A/" with (p, q) = j exist decouple from the system. Note, that 
the source term on the right hand side of (12.201) is known, since it consists of functions 
An,jm with either k = 1,. .. ,k — 2 or K = k — 1 and J7m G T^^i \ A/", which have been 
determined by the previous equations (I2.15P for 2, . . . , A; — 1. 

Then, since by (I2.20p now also the amplitudes ^fc-ij with j E M are determined, 
the right hand sides of the equations (I2.2ip are known, and we can determine by these 
equations the functions Afc.j™ for Jm E Tk \ Af, provided the system J\f of the pulses 
under consideration is closed under interactions up to order k, which guarantees Sj^ ^ 0. 
The fc-th order amplitudes A^j, j G A/", remain undetermined. Hence, performing the 
above procedure inductively up to A; = N, all equations (I2.15P or, equivalently, (I2.14p are 
satisfied, while no conditions are imposed on A^j. Thus, we obtain the following result, 
which establishes the formal derivation of the macroscopic dynamics. 

Theorem 2.1 Let {±{'dj,ujj) : j = 1, . . . ,1/} be a set of u E N different pulses, closed 
under interactions up to order N according to Definition \1.4\ and let the amplitudes A^j 



of the multiscale ansatz X^'^ given by (I1.13P satisfy the macroscopic equations (I2.17P for 
k = 1, (12.180 for k = 2, . . . , N—1, and set A]\rj = 0. Moreover, set the functions A^^j^, 
Jm E Tk \ Af, according to (12.200 for k = 2, (12.210 for k = 3, . . . ,N. Then, the ansatz 
X^'*^ constructed by these functions satisfies the microscopic model (II. 2p up to order , 
i.e. for the residual terms given by (127[3|) we have res(X^'') = 0{£^+^). 

However, all results obtained in this section are formal in the sense that they follow 
from the a priori assumption that solutions to (II. 2p retaining over time the form x = 
X^'*^ + 0{e^~^^) exist. Whether this is indeed the case, i.e. whether an approximation 
X^'*^ constructed by the solutions of the derived macroscopic equations stays close to an 
original solution of (11.20 over macroscopic time intervalls, is discussed in the next section. 
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3 Justification 

3.1 Estimate of the residuum 

In order to justify the macroscopic equations derived in the previous section, we will need 
an estimate of the residual terms res {X^''), see (EUS]), with respect to the f(T) -norm, 

11^11^2 = 1^71^- From fl2.13p we obtain 
76r 



N{N+1) 



|res(X;J'^)||,, <£^+it^+i||,2 + £^+lt;v+2||.2+ ^'\\"^\ 



k=N+l 



+ E 1 1 ( {^n') .^.-Kn I + 1 1 n+i {xin I 



provided of course the series over a G F exists. From the definition of tk in f l2.2p we obtain 

\\tN+iy + £\\tN+2\\i^ < Yl \\drAN,jJ\i2 + eo\\dlAN,jJ\e2 + \\dlAN-i,jJ\e2) 

(3.1) 



for £ G [0,£:o]) > (with A]^_i^j^ = when = 1 or m = A^). For given by (12.61) 
we first estimate (cf. (12. 4p ) 



\De,k-i,j„.,a\ < 2\ia-Vyf''Ae,jJ < 2\a\''-' Yl \DyAe,j^^ 

\a\=k-l 



(3.2) 



for £ = 1, . . . , A^, A; = £,... , A^+1 (with Ai^j^ = Af}^^ when k = N+1), where we used 



s d 



i=l di=l 
d 



di,...,ds=l i=l 



di ,...,d.,=l i=l i=l lcr|=s 



with a = (ai, . . . , a^) G and D^^ = d^yl ■ ■ ■ d^^, \a\ = (Ti+ . . . +ad. With ([S2D and 
we obtain 



<i+...+«n<fc-i «Gr i=l 



£2 



rt-1 



I I —kn~^n 



«l + ...+^„<fc i=l |o-|=A;i 
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with ||a;||£oo = max|x-y|, Ca,b,k,n '■= 2"^^ |a„^Q|(H-|a|'^ ") + and where we used that 
(n <)ii + . . . + in = k = ki + . . . + kn with £j < ki imphes ii = ki. Hence, by (12.61) we get 

N{N+1) N{N+1) N 

E ^'ik.ii,. E -0-"-' E E E E 

k = N + l k = N+l „-\ fc-l 1 I 1 fel + ... + fen = fe {Jm,,---,Jm„)eT^ li + ... + l„<k 

LjV + iJ-^-^ l<fcj<JV + l Jm-STi. m,<ei<k, 

n-1 

(n E E (3-3) 



I I — I cr I — 



for e G [0,eo] with Ca,b ■= ^I^ax^ (^2" ^ |a„^„|(2 + + and A^^j^, = A 



when fci = Finally, since by 14, We C^+^{R) 

\V^+iJx)\ < aN+i,a,.o\xf^\ \W'j,+,ix)\ < bN+i,.o\xf''' (3.4) 
for X G [— 2xo,2a;o], xq > 0, with aN+i,a,xoy bN+i,xo > (cf. (11 .4^ ). and since 

N N 

ll^^1U<^E^o-' E ii^'^.^Ji^-' \\^N%.<eY.^o-' E (3.5) 

fc=l JmeTfc fc = l JmeTk 

(cf. (11.13P ) for e G [0,eo]) we obtain, considering || (-^tv'^) .+q,||£2 = ||^7v''^||£2 ^^i a G F, 

E II^Af+l,a((^iv'^).+Q~^7v'^) IL2 + ||^iv+l(^^'^) 11^2 



< ^ aiv+i,Q,a;o 1 1 (X^'^) .^^-X^'^ 1 1^00 1 1 (-^^'^) .+Q~"-^^'^ 1 1^2 + bN+l,xo W^n'^ 1 1 £00 I |-^7V'^ 1 1^2 



TV ^ N 



k = l JrneTk k=l JmSTfe 



for e G [0,eo] and ^ < a^o with cv,w,xo ■= '2'^'^^^aN+i,a,xo + bN+i,xo- 

k=i Jmf^Tk aer 
Hence, for potentials 14, VF G C^+2(M) (a G F), cf. ([HD, which satisfy 



max ( 2, kn.al I'^r^ ) < ^ 

AT 

m, E 

<oo with Xo>J2^o E ll^'fc.^.JI^-'' (3-7) 
aer k=i JmeTk 

we obtain by (13. ip . (13. Sp and (13. 6p for Eq, tq > the estimate 

1 1 res (X;^-^) 11^, < e^+'Cr for e G (0,£o], t e [0,ro/£], (3.8) 
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with Cr > independent of e and t, provided the estimates 

\\9rAN,jJ\e2, Wd"^ An^ J J\e2, ||9^Aiv-i,j„.||f2 < oo, 

\\D^yAk,jJ\e^, \\D^yAk,jJU2 < oo foT\ar\=i-k,i = l,...,N + l,k = l,...,i (3.9) 

are satisfied uniformly in r = et e [0, To] (with Ak,j^ = Afl^^ when i = N + l,d. 1^). 

However, the k-th order amphtudes A^j, k = 1, . . . , N — 1, j = are obtained 

as solutions of the partial differential equations fl2.17p . fl2.20p . and the functions A^^j^, 
k = 2, . . . , N, Jm &Tk\ M, are given via the formulas (12.181) . (12.211) . All these functions 
depend on the (continuous) macroscopic time and space variables t = et and y = erj with 
1] = '-f for 7 G r. Thus, we have to relate the above i°°{r)- and £^(r)-norms to the norms 
of the solution spaces of the appearing partial differential equations. This is the purpose 
of the following lemma. 

Lemma 3.1 For d ^N, s > d/2 there exists a Cp > 0, such that 

||0(£(- + M)ll£2 <Cp£-'^/l0||H.(R.) 

for all (j) G ff(R'^), aeV, he[0, 1] and all e G (0, sq], Eq > 0. 

Proof: By Sobolev's imbedding theorem, there exists a > 0, such that ||v5||l°°(c) ^ 

c. ||<^||h»(C) for all ^ G W{C) with C := {higi + ... + hdQd ■ {hi, ...,hd) G [0, 1)^} (cf. 
(II. ip ). Hence, since ha = h{aigi + . . . + adgd) = [hai]gi + . . . + [had]gd + ^* with G Z, 
K G C, we get 

110 ie{- + ham% =J2\<P (^(7 + ha))f = |0 (.(7 + h^))f kJ^U (eil + -mUc) 
7er 7er 7er 

7Gr 

with Cp := Cs max{l, ^q}, obtaining the latter inequality by rescaling. □ 
Hence, under 

Assumption 3.2 For d E N and s > d/2 there exists a tq > such that the functions 
Ak,j,^ : [0,ro] x R'' ^ C, k = l,...,N, G %, given by (^J7^, (^M, <^M, (M) 
satisfy 

d, A^,j,^, dlA^,j^, 5^^jv-i,j,„, DlAu,j,^ G C ([0, ro]; W{W- C)) for \a\<N + l-k. 
we obtain 

Lemma 3.3 Under Assumption \3lM and for Eq > 0, Va,W G C^+^(]R) (a G T) as in 

N 

(ll.4p . (13.71) with xq > J2 ^0 ^Sl,x H^fc.j:^ ||oo, there exists a Cr > independent of e 

k=i JmeTk ^e[0,n)] 

and t, such that 

||res(XjJ'^)||^, < e''+^"'/^Cr forse (0,£o], t G [0,ro/e]. (3.10) 
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Proof: By Assumption I3.2[ Sobolev's imbedding theorem on M*^ with s > d/2, and 
Lemma [3.11 all norms in (13.91) are uniformly bounded for r G [0,ro], and (13.81) gives the 
assertion of the lemma. □ 



Remark 3.4 (a) Note that the order of e asserted in (13.101) could only be obtained due 
to the linear dependence of the estimates (13.31) . (13. 6p on the £^(r)-norms of DyA^^j^. 
Moreover, Assumption 13.21 implies by Sobolev's imbedding theorem the uniform bound- 
edness in time of these functions with respect to the L°°(]R'^)- and hence the £°°(r)-norm. 
Furthermore, assumption and theorem yield A^^j^^^r, ■) G C^"'"^~'^(]R'^), justifying opti- 
mally the Taylor expansion used in order to obtain (12.31) . (b) Using the bounds implied 
by Assumption 13.21 and the constants of Lemma [3TT] and of Sobolev's imbedding theorem 
on the constant Cr in fIXTUD could be given explicitly via dSH), (IS3D, dSSD, fIXTD . 

Remark 3.5 According to (12.181) . (I2.2ip and definition (12.101) of Sk-2,jm^ using the prop- 
erty A,B eC ([0,ro];ff(M^)) ^ AB e C ([0, tq]; (M^)) for s > d/2 (cf., e.g., jXF02l 
Th. 4.39]), it follows that Assumption 13.21 is satisfied iff 

Ak G ([0,ro]; (ff+l"l(M'^; C))'') for A + |(t| < + 1 - A;, A; = 1, . . . , iV - 1, (3.11) 

where := {A^^i, . . . , A^^u)'^ (cf. for A^ = 3 the examples in Sec. I4.1.3[ 14.3.31 below). 
Recall here, that A^f remains undetermined by the formal derivation procedure, and can 
thus be assumed as identically vanishing. 

The determining equations (I2.17P for Ai and (I2.20p for A^, k = 2, . . . , N — 1, are the 
(semilinear for k = 1, linear inhomogeneous for = 2, . . . , A^ — 1) symmetric hyperbolic 
systems 

d^Akir) + AkAkir) = /fe(r, A(r)) (3.12) 

where Ak '■ D{Ak) Xk is the infinitesimal generator of a Cq semigroup on Xk : = 
(^RM.-i^ ^d.QY ^i^j^ 3 /)(^^) 3 (H^*(R'^;C))'' =: n, Mk > l(cf., e.g., [RR931 
§11.3.1], |Gol851 §11.2.9]). Since fi does not depend explicitly on r and is quadratic, and 
hence locally Lipschitz continuous in Ai, for Mi > d/2 it holds /i : Yi — »• li, and we obtain 
by standard results of semigroup theory (cf., e.g., [Paz83l Th. 6.1.7], |Gol85t §11.1.3]) that 
there exists a Tmax < oo such that (I3.12p has for /c = 1 and Ai{0, ■) G Yi a unique classical 
solution Ai G C^([0, tq]; Xi) n C([0, tq]; Yi) for tq < r^ax- For A; = 2, . . . , A^ - 1 we have 

fk{T, Akir)) = F{r)Ak{r) + Gir)Mr) + 9k{r) 

with F,Ge C([0,ro]; {}1^^'{R'^;C)Y'"') and gk e C ([0,ro]; (H^^^-i-2(M'^; C))'') . In order 
to obtain unique classical solutions Ak G C^([0, tq]; Xk) H C([0, tq]; Yk) of the initial- value 
problems (13.120 with ^4^(0,-) G Yk, we need that fk : [0,Tq] x Yk Yk is continuous 
in (r, Ak) and uniformly Lipschitz in Ak- Hence, we need Mk-i > Mk + 2, and thus 
Mk > Mn-1 + 2{N-k-l) for A; = 1, . . . , A^-1. Moreover, according to fl3?TT|) . we need 
Mfc>s + A^ + 1-A;, and in particular M^-i > s + 2. Hence, for Mk = s + 2{N-k) 
all conditions on Mk are satisfied, which means that assuming initial data ^4^(0, ■) G 
ff+2{JV-fc)(M'i;C), for the equations fl^TTD and ^TM) we can guarantee flXTTD . and thus 
Assumption 13.21 

In particular it is necessary and sufficient to assume Ai(0, ■) G H'^"*"^*^^~-'^^(R'^; C). Then, 
we can assume ^1^(0, ■) = for A; = 2, . . . , A^ — 1. This does not yield identically vanishing 
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k-th order amplitudes A^, due to the source terms in fl2.20p . Of course, also some of 
the first-order amplitudes can be chosen as initially vanishing. For example, setting in 
the case of the three-wave interaction fl4.14l) y4i 2(0, ■) = 0, and assuming that Ai i{0,-) 
and ^1^3(0, ■) have disjoint supports, the amplitudes Ai^i and Ai^^ will be transported by 
their different group velocities V,jn('(9i) and V^n(t93). The moment they interact, i.e., 
when their supports intersect, the amplitude 2 arises. Note, however that if two of the 
three amplitudes are assumed to vanish initially, then the third one is just transported, 
and does not give rise to other amplitudes, since it does not interact with other pulses. 
(Confer also the discussion on generation of pulses in Sec. Hlbelow, in particular Sec. 14.3.11 
Case 121 and Sec. liXTl ) 



3.2 The justification result 



Having obtained the estimate fl3.10p of the residuum, we are now able to establish the 
main result of our work, namely the justification of the macroscopic equations (12.171) . 
(I2.18p . (I2.20p . (12.211) obtained by formal derivation in Section [21 as giving the effective 
dynamics for an arbirtrary number of amplitude-modulated pulses in multidimensional 
lattices with scalar displacement. More precisely we show: 



Theorem 3.6 Let T be the d- dimensional lattice f lLTjl (den) and let V^, W E C^+'^( 
(a E T, N E N, N > 1+^) be the interaction and on-site potentials (11.41) of the micro- 
scopic model (II. 2p satisfying Va{x) = 1/_q,(— x), (II. 8p . (13. 7p . Let {±{{}j, Uj) : j = 1, . . . , z/} 
be a set of u E ^ different pulses, closed under interactions up to order N according to 
Definition \1.4\ and let the functions A^^j^ : [0,tq] C of the approximation X^'"^ , (11.130 . 
solve the macroscopic equations (12.170 . (12.180 . (I2.20p . (I2.2ip as described in Theorem \2.1[ 
and satisfy Assumption \3.B for some tq > 0. 

Then, for each c > there exist eo,C > such that for all e E (0, £0) one? any solution 
X of ([L2D with 



(x(0),x(0)) 



(X^f,(0),X;J'f,(0)) 



(3e{i,n-I], 



(3.13) 



it holds 



{x{t),x{t))-{xi'l,it),Xi'l,it)) <Ce^ for tE[0,To/6]. (3.14) 

Remark 3.7 In view of Remark 13.51 the Assumption 13.21 above can be replaced by the 
assumption AkjiO, ■) E H^+2(^-^)(M'^; C) for A; = 1, . . . , _ 1, j = 1, . . . , zy with s > d/2 
(and in particular by Aij{0, ■) E H"+2(7v-i)(^]gd. Akj{0, ■) = 0, k > 2), which moreover 
guarantees the solvability of the respective amplitude equations (I2.17P for = 1 and 
fl2:20|) for > 2 up to some tq > 0. 

Proof: We write the microscopic model (II. 2p as a first order system in Y := ^^(r)x£^(r). 



X = Cx + A4(x) with X :- 



X 
C 



Mix) :-- 





M{x) 



(3.15) 
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where J : ^^(r) ^ ^{T) is the identity, and C,M : ^^(r) ^ i'^{T) are given by 

{Cx)y := ^ai,a(x^+„-x^) := ^ V^a^^lx^+a-x^,) - W^ix^) (3.16) 

with ai^a, y^a^ ^1' ^2 in fll.4p . On the Banach space Y we use the energy norm 

^l,a I i2 , I I |2 



defined in such a way that its square is twice the harmonic part of the Hamiltonian 7i, 
cf. (11.61) . As is well known, the fiow of the linearized system x = Cx preserves this norm, 
i.e. its associated semigroup e*^ satisfies ||e*^||y^y = 1 (cf., e.g., |GM04l Prop. 3.1]). 
Moreover, since || ■ \\i2 and || ■ H^; are equivalent by the stability assumption (II. Sp . 

/^-||a;||£2 < < /i+||x||£2, < /i_ := min 0(79), := maxfl^-d), 

which follows by Fourier transformation, we obtain the equivalence of || ■ \\fiy^£2 and || ■ ||y, 

/I_||x||^2x^2 < \\x\\y < Jl+\\x\\e2^(,2, Ji_ := min{/i_, 1}, /!+ := max{/i+, 1}. (3.17) 

We show that the error Re{t) := e~^{x{t) — X^'^(t)) between a solution x of (13.150 
and the approximation X^'^ := satisfies 



||^e(0)||y <Co =^ \\R,{t)\\y<D for £ < eo, £t < To. (3.18) 

By (13.171) . this is the assertion of the theorem with Xfj'^ instead of X'^'ti- However, 
according to (11.130 . Lemma [3.11 and Assumption 13. 2[ there exists a cjv > such that 

1 

\\X^^\t)~X^tAt)\\Y < E + {4^rAM,j4,, + \^jj\\AN,j4,,)y 

< e^-'^'^cn for e < Eq, et < tq. 
Hence, by ( KT7} and /3 < iV - d/2, we obtain from fl3:T3D 

p.(0)||^<£-'^||2?(0)-X;}f,(0)||^ + £-^||X^'^(0)-Xi5'!i(0)||^ 

^ ~ , N-d/2-l3 r ^ 

< fl^c + Eq Cn =: Co for e < Eq. 

Vice versa, from the right hand side of (I3.18P and f3 < N — d/2 we obtain (13.140 

E-q^t)-Xi'l,it)l,^^,<Jlz'(^\\RM^^ 

< JlZ^^D + Eq~'^^'^~^cn) =■ C for e < Eq, Et < tq. 

It remains to prove (13.181) . Inserting into (13.151) we obtain the differential equation 

R, = CR, + E-^ [m{X^'' + e^Re) - M{X^^') - res(X^'^) 
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with res(X^'^) := X^'"^ — CX^'^ — and write it in its integral form 

e*-^ the semigroup to = CR^. Since ||e*'''||y^y = 1 (cf. above), taking the norm gives 
\\Rsit)\\y<\\RM\\y + e-^ l\\\M{xi'\s)+e^R^^^^ 

+ ||res(X^'^(s))||^)ds. (3.19) 

We estimate the norms on the right hand side of (13.191) . According to the formal 
derivation in Section [2], X^'"^ — CX^''^ — Ai[X^''^) = res(X^''^). Hence, Lemma 1331 gives 

||res(X^'^)||^ = ||res(X^'^)||,, < e^+i-'^/^C, for e G (0,£o], t e [0,To/e]. (3.20) 
Next, we estimate the norm of the nonlinear terms. From (I3.16P we obtain 

||A^(x) - M{y)l, < J2 \\VUx.+^-x) - Vl^iy.+^-y)l, + WW^ix) - W^iy)^, 

for X, y G ^^(r), where we denote {x.+a)^ '■= 2^7+0 and (y{x))^ := V{x^) for V : R ^ W. 
Since 14, W G C^+2(M) with X > 2 and V^'M = 0{\x\), W^'{x) = 0{\x\) (cf. ([L4])), we 
obtain by the mean value theorem 

KM-VlM\<c^,,,{\x\+\y\)\x-yl \W^ix)-W^iy)\<c,,{\x\ + \y\)\x-y\ 

for X, y G [— 2xo,2xo], Xq > 0, with Ca,xo:(^xo > 0- Hence, there exists a cm,xo > (cf. also 
(13. 7p ). such that for x,y & Y 

||A^(x) - A^(y)||^ < CA^,a;o (lklk°° + l|y|k->) p- y|k for llxll^oo, \\y\\eoo<xo. (3.21) 

Now, we apply (I3.2ip on x : = X^'' + e^R^ and ^ := X^'". By ([33]) and Assumption [3j 
there exists a Cyi > such that 

||X^'^||^^ < ioi £ < £t < Tq. 

We set D := (cq + ej'"'^/^"^ tq a)e^o^^ with co,C,. as in fl3J8|) . fICTD . Cm ■= ^cm,xoCa, 
and £0 > such that £oCa < Xo/2 and £Q^^D/fi_ < ca, thereby assuming (3 > 1. Since 
||-Re(0) lly < Cq < D and ||i?£(t) ||y is continuous, there exists for every £ G (0, ^o] a tf^ > 0, 
such that ||-Re(t)||^ < D for t < t'j^. Then, it holds 

ll^iv'ILoc, WX^'' + ^'^RsW^oo < Xo (OT £ < £0, t <mm{To/£,t'jj}, 

and (I3.2ip gives 

\\M{xi''+£^R,) - M{xi'')\\Y < £^+'CM\\Re\\Y for ^ < ^0, t < mm{To/£,t%}. 
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Inserting this estimate, fl3.20p and ||-R£(0)||^ < Cq into (13.191) . we obtain 

||-R£(^)||y< (co + Eq''^^'^''^ ToCr^ +eCM ||-Re(s)||yds ioT e < Eo, t <mm{To/e,t%}, 

and Gronwall's lemma gives 

\\Re{t)\\y< (co + To a) e^*^^ foTe<6o,t<mm{ro/6,t%}. 

Hence, < D ior e < eq, t < tq/e < tf^,, and the proof of (13.181) . and thus of the 

theorem, is completed. □ 

4 Examples 

In the present section we give, on the one hand, a complete classification of all possible 
types of macroscopic systems for the first-order amplitudes Aij of z/ = 1, 2, 3 modulated 
different pulses, which form sets {±('t?j, tu-,) : j = 1, . . . , z/} closed under interactions up to 
order N = 2. The classification for a given number of pulses results from the number and 
type of resonances (of order 2) between them. These resonances can be traced back on 
the coupling terms of the first-order amphtude equations. Each case leads to a different 
macroscopic system, which can not be obtained from another even when the resonance 
conditions for the former are a special case of those for the latter (cf. the discussion in 
Sec. 14.3.11 CaseH]). In particular, additional resonances lead to additional terms in the 
macroscopic systems (cf., e.g., in Sec. 14.3.1] the Cases [3] with H] or [21 O and [2] with that 
of Sec. I4.3.2p . Naturally, the more pulses we consider the more different constellations of 
resonances can appear, which of course can include also closed subsystems (cf. e.g., (14.21) 
with (112!) or Sec. lOAl Case[Il and with Sec. lOm Case [3]). More precisely, withm 
a system of v pulses with < uoi < . . . < uj^, the resonances of order 2 

{'dj+^i,ujj + uji) = {'dk,i^k), k>j,i, k,j,i = 1,. . . ,u, 

can appear. 

On the other hand, we exemplify on the considered sets the essence of the closedness 
condition. Only when it is satisfied we are able to obtain complete effective dynamics, 
in the sense that only then we can choose (except for regularity restrictions) arbitrary 
initial data for the macroscopic amplitudes of the pulses (cf. in particular the discussions 
of non-2-closed sets in Sec. 14.1.21 and 14.2.21 and of non-3-closed ones in Sec. I4.1.3p . 

Recall that whenever in the following a set is closed up to interactions of order N = 2 
(2-closed) the derived macroscopic equations provide us with the first-order amplitudes 
Ai j and the functions A2^j^, Jm G 7^ \ A/" of the second-order approximation X^'^ such 
that res(X^'^) = 0{e^\ while the second-order amplitudes A2^i remain undetermined, 
enabling us to choose A2^j = 0. Thus, according to Theorem 13.61 and Remark 13. 7[ for 
one-dimensional lattices (oscillator chains) with Va,W & C^(]R) the given macroscopic 
equations for Aij with Aij{0,-) G II^(]R) establish the effective dynamics of the corre- 
sponding modulated pulses in the sense of Theorem l3.6t i.e. with A^— 1 = 1 and (3 G (1, 3/2] 
in (I3.13P and (I3.14p and for macroscopic time intervalls [0,ro], where tq > is limited 
only by the time of existence of the macroscopic solutions Aij. 
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For deriving valid effective dynamics in multidimensional lattices with d > 2, however, 
we need the considered set of pulses to be closed under interactions up to order > 
1 + d/2. In particular for d = 2 or 3 we need = 3. Then, for Va,W G C^iR), the 
equations determining the macroscopic coefficients of X^'' with Aij{0, ■) e R^iW^) and 
^2^(0, ■) G H'*(M'^) (or y42j(0, ■) = 0) give the effective dynamics of the corresponding 
modulated pulses in the sense of Theorem 13.61 i.e. with — 1 = 2 and j3 G (1,2] or 
P G (1, 3/2], respectively, in fl3.13p and (13.141) . We give in the following the corresponding 
macroscopic equations in the case of a single pulse (Sec. 14.1.1^ Sec. I4.1.3P and in the case 
of a three-wave-interaction (Sec. 14.3.11 Case [2], Sec. I4.3.3p . 

We would like to remark that of course each single model presented in the following 
could be discussed and interpreted more extensively than done here. Note, in this context 
that the presented macroscopic equations are (at first formally) valid for pair-interaction 
potentials among atoms at an arbitrary distance in lattices of arbitrary dimension. Spec- 
ifying the potentials (as for instance by considering only nearest-neighbour interactions) 
and the lattice dimension would lead naturally to more concrete results and possible inter- 
pretations. Nevertheless, we believe that the comparative presentation of all these models 
is helpfuU in getting an impresion of the variety of possible interactions of pulses, and 
even may serve as a reference for further investigations. 

Concerning the following examples the most crucial question is of course whether 
the resonance conditions corresponding to each model can be satisfied at all for a given 
potential in a given dimension. Here, we restrict ourselves in addresing this question 
exemplarily in Sec. l4.4l for the case of a three- wave-interaction in an oscillator chain [d = 1) 
with nearest-neighbour interaction potential and a stabilizing on-site potential. This 
example shows that the satisfaction of resonance conditions is equivalent to the solution 
of algebraic equations in n variables, where n depends on the number of independent 
pulses, cf. also |GHM08al Sec. lll.E]. In the present case n = 2, and it turns out that a 
three-wave-interaction can appear only when the leading, harmonic part of the interaction 
potential is repulsive. Recall here, that the satisfaction of resonance and nonresonance 
conditions depends only on the harmonic parts of the involved potentials, since they 
are linked to the dispersion function of the (linearized) lattice, see (11. 7p and Definitions 
11.31 and II. 4[ Moreover, note that, since closedness (nonresonance) conditions appear as 
avoided solutions of algebraic equations, as shown in our example, they are usually more 
easily satisfied than resonance conditions. However, aside from the observation that the 
latter depend on the number of involved pulses, the range of the interaction potential and 
the dimension of the lattice, we do not follow here these questions further. 

4.1 A single pulse 

We start with the case z/ = 1 of a single pulse ±{i}i,uJi) with uJi = > and consider 

the second-order approximation fll.l3p 

X^^' = £Ai,iEi + (v42,iEi + A2,(i,i)E2 + iA2,(i,_i)) + c.c. (4.1) 

Here, E±i(t,7) = eK'^±it+^±i-y) with u;±i = ±uji, ^9±i = ±^9i. The muhi-index (1,-1) G 
A^^ is the representant of ±(1, —1). The factor i in front of ^2,(1,^1) G M arises, since we 
abbreviate the complex conjugates of all explicitly written functions on the right hand 
side of (14. ip by c.c. 
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4.1.1 2-closed system of a single pulse 

According to Definition II ■4[ the set of pulses {±('(9i, tui)} is closed under interaction up 
to order 2 if 

the pair (2'(9i, 2ti;i) does not characterize a pulse. 

Recall here, that the pair (0, 0) corresponding to the functions E-|-(i _i) = Ej-iEzpi = 1 
does not characterize a pulse by our stability assumption (11.81) : > for all -i? G Tp. 

In particular, (-dijUJi) ^ (0,0), which yields (2'i9i, 2c<Ji) ^ ±('(9i,u;i), meaning that the 
functions E^j^ generated by the pulses E-ti can not equal one of the latter. Hence, here 
the closedness condition is equivalent to the nonresonance condition 

5(1,1) =fi2(2i9i)-4a;iV0, 
and we can determine the function A^^x in (14.11) by the equation (I2.17P 



9.Ai,i-V^l](?9i)-V,Ai,i=0, (4.2) 
and the functions A'i^{x,±\) by the equations (I2.18P 

/l2,(l,l) = T /ll,l, ^2,(1,-1) - |/li,i| 

(^(1,1) Ol 



with 



C(i,i) = -4i ^ a2,a sin^ ( ) (^i'*^) ~ ^2- (4.3) 



According to (14.21) . the amplitude A^^i of a single modulated pulse, which does not 
generate a further pulse via self-interaction, is simply transported with its group velocity 
V^fi('i9i). Note here, that this does not mean that the pulse ±('(9i,ci;i) could not possibly 
interact with some other pulse ±(-(92; "^2)- It just means that if initially all other pulses 
have vanishing amplitude, they cannot be generated (i.e., exhihit at a later time a non- 
vanishing amplitude) by a single pulse, except in the case of self-interaction, = 0, 
which is discussed in the following. 



4.1.2 Non-2-closed system of a single pulse 

The closedness condition for the system of a single pulse is violated if 5(i,i) = 0, i.e., in the 
case where the pulse ±{{}i,uji) generates via self-interaction a further pulse ±{2{}i,2u!i). 
Considering still, by the ansatz (14. ip . only the pulse ±{i)i,uJi) we would obtain Ai^i = 
(cf. the equation for v42,(i,i) in Sec. 14.1.11 above). This means that we would not consider 
any modulated pulses at all, and in particular that non- vanishing initial data y4i,i(0, ■) 7^ 
would be excluded. The situation is resolved if we take into account also (at least) the 
generated pulse {^2,^2) '■= (2'(9i, 2ti;i). Since {^2,^^2) 7^ ^(^i,'^!), this leads to the case of 
a set of 1/ > 2 pulses, which for z/ = 2, 3 is discussed in Sections 14.21 and |4^ respectively. 
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4.1.3 3-closed system of a single pulse 

The third-order approximation fll.lSp for a single pulse ±{'di,ui) with ui = > is 
given by 

X^'' = X^'' + 6^ (As.iEi + A3,(i,i)E? + A3,(i,i,i)E? + ^A3,(i,_i) + c.c.) (4.4) 

where X^'^ is the second-order approximation (14.11) . Here, the set of representants is 
% = {±1, ±(1, 1), ±(1, 1, 1), (1, -1)} with T2 = {±(1, 1), (1, -1)} and Af = {±1}. Hence, 
J\f is closed under interactions up to order 3 if the nonresonance conditions 

5(1,1) = ^^'(2^i) - 4ujI ^ 0, 5(1,1,1) = ^^'(3^i) - 9ul ^ 

are satisfied. The second-order amplitude ^2,1 is determined by equation (12.201) for A; = 3 

5.^2,1 - V^fi(^i)-V,A2,i = — ^ I (2^i^ + 4^ + 3c(i,i,_i)) |Ai,i|Mi,i 

2ifi(?9i) \ V 5(1,1) &i I . , V ' ' ' 



+ ^ E «i,ae^"-"(«-V,)Mi,i - 9^Ai,i ) . 



5(1,1)^3,(1,1) = - 4ia;i5, + 2in(2t9i)V^r](2^i)-Vj,) ^2 , + 2c(i,i)Ai,iA2,i 

"(1,1) 



and equation (I2.2ip for A; = 3 gives 

&i^3,(i,-i) = -262^1,1^ + 2Ai,i ^ a2,„ (l-e-'''i-°) a-Vj,Al7 + c.c, 

C(i,i) 
<^(i,i) 

+ 2^a2,„ {f^^ '--^^-^) a.Ai,iV,Ai,i, 

/ C(i 1) \ 3 

5(1,1,1)^43,(1,1,1) = ( 2c(i,(i,i))— ^ + C(i,i,i) ]A^ ^, 
^ ^(1,1) ^ 

where we used the functions v42,(i,i) and ^2,(1,-1) calculated in Sec. 14.1.11 and where C(i,i), 
C(i,(i,i)), and C(i,i,i), C(i,i,_i) are given by (14.31) . (12.81) . and (12. 7p . respectively. Hence, 
the second-order amplitude ^2,1 is determined by an inhomogeneous transport equation, 
where the source term is determined by the first-order amplitude Ai i calculated by (14.21) . 
Then, since 5(i,_i), 5(i,i), 5(i,i,i) 7^ by the nonresonance conditions, and knowing the 
amplitudes Ai i, ^2,1, we can calculate the functions v43,(i,±i), ^3, (1,1,1). The third-order 
amplitude ^3,1 remains undetermined. Setting ^3,1 = 0, and inserting the above functions 
and X^'^ from Sec. HXD into (|13D, we obtain X^'' with res(X3'^'") = O(e^). 

Note here, that this derivation of amplitude equations for ^2,1 and A3,(i,±i) is possible 
only since we anticipated for functions Aa,.,j-„ with J7m G ^ and not only with G 
for /c = 2, 3 in (11.131) . In the latter case, we would have assumed ^2,1 = ^3,(i,±i) = 
identically, and would obtain instead of the above equations for ^2,1 and ^3,(1,^1) addi- 
tional conditions on the fist-order amplitude 74i,i which would at least restrict our choice 
of initial conditions. 



26 



The 3-closedness of the system would be violated if (^(1,1,1) = 0, i.e., if (S"!?!, Scui) is a 
pulse (resonance of order 3), not considered by our system of a single pulse. Ignoring this 
generated pulse would give 1 = identically (cf. the equation for ^3,(1, 1,1)), prohibiting 
non-trivial dynamics even for the first-order amplitude Ai^i. Again, the remedy consists 
in taking into account also the generated pulse, thus considering a set of two pulses. Note, 
however, that the resonance = is detected only by the third-order approximation 

X^'"^. Thus, even in its presence the system remains 2-closed (when 7^ 0), and the 
derivation of valid effective dynamics of the first-order amplitude Ai^i is still possible for 
one-dimensional lattices. This is in line with the order f3 G (1,3/2] of the error in (13.131) 
and (13.141) . which for = 2 is too coarse to detect this interaction. Of course, if one is 
interested in valid higher-order effective dynamics (with > 3) the same problems as 
discussed above arise. 

4.2 Two pulses 

We consider the case z/ = 2 of two different pulses {■di,uji) 7^ (^92, ^^2) with Uj = > 0, 

j = 1,2. The second-order approximation (I1.13P reads 

2 2 

X^^' = e J2 ^i,.E, +e\Yl (^2jE, + A2,u,)E^) + A2,(i,2)EiE2 + A2,(i,-2)EiE_2 

+ ^^,(1,-1)) +C.C., (4.5) 

with Ej{t, 7) = e'(^^*+''^-^), j = ±1, ±2, where uj_j = -ujj, = -dp Here, (1, 2), (1, -2) 
represent also (2, 1), (—2, 1), respectively, and (1, —1) represents ±(1, —1), ±(2, —2). 

4.2.1 2-closed systems of two pulses 

According to Definition 11.41 the set of pulses {±('i9j, cUj) : j = 1,2} is closed under inter- 
action up to order 2 if 

the pairs {2'dj,2ujj), j = 1,2, (^9i±'(92, 1^1^1-1^2) do not characterize pulses, 
except if either (2'i9i, 2ci;i) = (792, (^2) or (2792, 2072) = {{}i,u!i). 

The presence or not of one of these two exceptional cases, leads to two different cases 
among 2-closed sets of two pulses: 

Case 1: Non-interacting pulses. 

Excluding the exceptional cases above, the closedness condition is equivalent to the non- 
resonance conditions 

= ^2(2^^.) - Alu] ^0, j = 1, 2, 5(i,±2) = ^^idi±d2) - (^i±^2)' 7^ 0, (4.6) 

and we can determine the first-order amplitudes Aij in (14.51) by the equations (12.171) 

v^n(79i)-v,A,i = 0, 

(4.7) 

5.Ai,2 - VM^2)-S/yAi^2 = 0, 
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and the functions ^2,0^) (j = 1,2), A2,(i,±2), ^2,(1,-1) by f!2.18p 



"^2,(jj) — 7 ^17) ^2,(1,±2) — -J ^1.1^1,±2, ^2,(1,-1) — ^^ U^l.ll + l^l,2| j 

with C(p_q) given by (12.191) . 

Note, that here all microscopic patterns within the e^-term of (14. 5 p are mutually 
different, except when either (3'(9i, Scui) = {'&2,^2) or {{}i,ui) = (3'i92, 3u;2), where we 
obtain only one respresentant for either {(1, 1), (—1,2)} or {(2,2), (1, —2)}, respectively, 
and the corresponding amplitudes are given by either ^2,(1,1) + ^2,(1,-2) or ^42, (2, 2) +^2,(1,-2) 
with the values from above. 

Hence, in the present case the two pulses do not interact with each other in leading 
order. Their first-order amplitudes pass through each other travelling with the respective 
group- velocities of the pulses they modulate, i.e., the corresponding transport equations 
are uncoupled. Thus, ignoring one of the two amplitudes (i.e. assuming it vanishes iden- 
tically in space and time) the other remains unchanged, and X^'^ is the same as for a 
single 2-closed pulse, see Sec. 14.1.11 



Case 2: Interacting pulses. 

In the case (2'i9i, 2c<Ji) = {1^2, UJ2) we have = 62 = ^(1,-2) = ^-i = 0, the closedness 
condition (14. 6 p reduces to the nonresonance conditions 

(5(2,2) = - 16ul ^ 0, 5(1,2) = ^^'(3^i) - 9ljI ^ 0, (4.8) 

and the equations (12.170 read 



drA,,2 - V,fi(2^i)-V,Ai,2 = -i^^^li 



(4.9) 



with C(i,i) given by 1^. Since T2 \ TV = {±(2, 2), ±(1, 2), (0, 0)}, we obtain by (l2A8l) 
the functions ^2,(2,2), ^2,(1,2), ^2,(1,-1) as in Case 1, whereas the functions ^2,(1,1), ^2,(1,-2) 
are included in ^2,2, ^2,1, respectively, which remain undetermined. (Compare also the 
values of the former functions in Case 1 with the right-hand sides of (14. 9p .) 

Here, the evolution equations for the amplitudes are coupled, in an asymmetric way 
that reflects that the pulse (^92,1^2) is generated by self-interaction of {{}i,u!i). In partic- 
ular, setting identically Ai i = 0, (14.90 becomes a homogeneous transport equation for 
Ai^2, whereas, setting Ai^2 = 0, (14.90 yields Ai^i = 0. The latter is exactly the situation 
in a non-2-closed system of a single pulse, see Sec. 14.1.21 In the contrary, we can describe 
the dynamics of such a self-interacting pulse by ( 14.90 . by taking into account via the 
ansatz (14. 5 p also the amplitude Ai^2 of the generated pulse, provided of course the system 
is 2-closed. (Note in particular how even for ^41,2(0, ■) = a non-zero Ai,2 is generated 
by (14.90 . when Ai,i(0, •) 7^ 0.) The situation of non-2-closed systems of two pulses is 
discussed in the following. 
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4.2.2 Non-2-closed systems of two pulses 
Case 1: Non-interacting pulses. 

The closedness condition fl4.6l) is violated if one of the pulses generates via self-interaction 
a third, {{^^jU^) = {2{}j,2u!j), j = 1,2, or if the two pulses interact to generate a third, 
(t^a, cus) = (-i?! ± ^2, "^1 ± ^^2)- In the present case neither of these pulses equals ±(^9j, ujj), 
j = 1,2. Of course, more than one of the four (or three, if {3'dj,3uj) = {{}i,uji), j 7^ i) 
corresponding nonresonance conditions ( 14. 6p can be violated, and moreover, even if only 
one new pulse arises it is not at all clear that the three pulses form a 2-closed set. 

Hence, we restrict our discussion to two typical cases where only one of the four 
conditions (14. 6 p is violated. In the first case, we assume 5(2,2) = 0. Then according to the 
equation for ^2,2 in Sec. 14.2. H Case 1, it follows Ai^2 = (if C(2,2) 7^ 0, e.g., when 62 7^ 
in (I2.19P ). Then, as discussed there (and setting ^2,(2,2) = 0), we can still describe the 
dynamics of Ai^2 as in Sec. I4.1.1[ However, we are not able to determine any non-trivial 
amplitude A12 for the self-interacting pulse ±{'&2,uj2) (for instance by assuming initially 
^1,2(0,-) 7^ 0). A first step to remedy this is to include the pulse {{l^yiVs) = (2'(92, 2ti;2) 
into the set of considered pulses, see Sec. 14.3. H Case [31 

In the second case we assume 5(1,2) = 0, which implies Ai^iAi^2 = 0, see the equation 
for ^2,(1, ±2) in Sec. 14.2.11 C. 1. Hence, Ai^i and Ai^2 should be for all times disjointly 
supported. Except for Ai^i = or Ai^2 = identically, leading again to the case of a 
single pulse, this violates the essence of interaction of the corresponding pulses, which 
thus can not be explained. (Apart from restricting the choice of allowed initial data for 
the amplitudes, this fails to predict their evolution e.g. if they move towards each other.) 
The correct set-up is once again to take into account also the generated pulse, see Sec. 
l4Xn Case [3 

Case 2: Interacting pulses. 

The closedness condition ( 14.81) is violated if 5(2,2) = or 5(i,2) = 0, generating a third 
pulse, {'&3,LJ3) = (4i9i, 4a;i) or (■(?3, cus) = (3'(9i, 3a;i), different from ±(-(9^, cjj), j = 1,2. The 
respective 2-closed systems of three pulses are discussed in Sec. 14.3. H Cases HI [51 Here, 
^(2,2) = yields as above Ai,2 = 0, but moreover, Ai^i = by (14.90 . Thus, noting that 
(1^2, '^2) = (2191, 2ci;i) is generated via self-interaction, we observe the same phenomenon 
mentioned already in Sec. 14.1.21 and Sec. 14.2.11 Case 2: 

if we ignore a generated pulse (by assuming it has an identically vanishing amplitude), 
then we have to ignore also (at least) one of its generators. 

In the second case, 5(i,2) = 0, we obtain Ai^iAi^2 = as in Case 1. Apart from the 
implications discussed there, here the macroscopic system (14. 9 p shows clearly an inherent 
direction in the generation of pulses: while Ai^2 = yields Ai^i = 0, assuming Ai^i = 
we still can describe the dynamics of Ai,2. Thus, 

generating pulses can be ignored, while generated ones can not. 

Hence, in order to obtain relevant macroscopic dynamics for a system of modulated pulses 
we have to take into account also the pulses generated by them, i.e., to require that the 
set is closed with respect to interactions, motivating Definition II. 4[ 
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Concerning the present case, when two pulses generate a third, new one, the appro- 
priate setting for describing their dynamics is given by considering the set of the three 
pulses, provided of course it is closed. We do this in the following. 

4.3 Three pulses 

We consider the case u = 3 of three different pulses {{}j,Uj) ^ {{}i,u!i) for i ^ j and 
2,j = 1,2,3, with ojj = > 0. The second-order approximation fll.l3p is given by 

3 3 

X^'" = £^1 ( 5Z (^2,iEj + v42,(jj)Ej^) + ^ (v42,(j,i)EjEi + A2,(j _»)EjE_i) 

j=i j=i i<i<«<3 

+ ^^(1,-1)) +C.C., (4.10) 

with Ej{t,-f) = e'(^^*+''^-T), j = ±1,±2,±3, where u.j = -uj, = -^j. Here, (1,±2), 
(1,±3) and (2, ±3) represent also (±2,1), (±3,1) and (±3,2), respectively, and (1,-1) 
represents ±(j, — j), j = 1,2,3. Else, we start from the assumption that {j, (j,j) : j = 
1, 2, 3} U {(1, 2), (1, -2), (1, 3), (1, -3), (2, 3), (2, -3)} is a set of (different) representants. 

In Section 14.3.11 we give the explicit evolution equations for the first-order amplitudes 
Aij, j = 1, 2, 3, and the functions ^2,^™, J7m G 72 \ A/" of the approximation (I4.10p for a 
set which is closed under interactions up to order N = 2. The equations depend on the 
presence, the number and the kind of interactions within the set. 

In Section l4.3.2l we consider a non-2-closed set for a three- wave-interaction and give the 
equations for the first-order amplitudes of a 2-closed system of four pulses with two three- 
wave-interactions. Finally, in Section 14.3.31 we determine a third-order approximation 
X^'^ for the 3-closed system of the three- wave-interaction of Sec. 14.3.21 Case [2l 

4.3.1 2-closed systems of three pulses 

According to Definition 11.41 the set of three different pulses {±('i9j, luj) : j = 1,2,3} is 
closed under interactions up to order 2 if 

the pairs {2'dj,2ujj), {'dj±'di,ujj±u!i) , j ^ i, (j,i = 1,2,3) do not characterize pulses, 
except if {^j+^i, ujj+uJi) = {§k,^k), k j,i. 

Since the latter exceptions correspond to interactions among the three pulses, we can clas- 
sify the possible settings according to the number of interactions and the pulses involved, 
up to relabeling, as follows: 

1. No interactions. 

2. Only a three-wave-interaction: (-i^i ± '62, oJi ± UJ2) = ("i^a, c^s). 

3. Only one self-interaction: (2^91, 2ct;i) = {1^2, ^2)- 

4. Two self-interactions: (2^91, 2u;i) = {'&2,uJ2), (2^92,2^^2) = (4'(9i,4cji) = {'&3,uj3). 

5. One self-interaction and one three-wave-interaction: (2'i9i, 2c(Ji) = (-(92,1^2), 
i^i + ^2,uJi + uj2) = (3^9i,3a;i) = (^93,^3)- 
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We assume that the nonresonance conditions given in the following are independent. 
It is possible that for special choices of pulses they can be reduced, cf. the example of 
Sec. 14.2. H Case 1. However, this does not affect the first-order amplitude equations, but 
only the coefficients of 'non-pulses'. Hence, the following classification remains unaffected. 

Case [It No interactions. 

The closedness condition is equivalent to the nonresonance conditions 

%i) 7^ 0' j = 1'2,3, 5(i,±2) 7^ 0, 5(i,±3) 7^ 0, 5(2,±3) 7^ 0, 
with 5(j-±i) = n'^{^j±^i) - {ujj±ujiy for j,i = 1,2,3, and OTTl) yields 

'a,Ai,i-V^fi(7?i)-V^Ai,i = 0, 

< 9,Ai,2- v^n(i?2)-v^Ai,2 = 0, 

By (EH]) we obtain (j, z = 1, 2, 3) 

Since the pulses do not interact, the amplitudes Aj are transported independently of each 
other by the respective group velocities 'W^VL{dj). If their trajectories intersect, they sim- 
ply pass through each other. Since the equations are uncoupled, neglecting one or two 
amplitudes by setting identically Aij = 0, yields the corresponding systems for two or 
one pulse, see Sec. 14.2. H Case 1, and Sec. 14.1.11 

Case Et Only a three- wave-interaction. 

We assume that the resonance condition 

(t9i+?92,^i+^2) = (^93,^3) (4.11) 

holds, and that self- interactions are absent. The closedness condition is equivalent to the 
nonresonance conditions 

%,)^0, J = 1,2,3, 5(1,3)^0, 5(i,_2)^0, 5(2,3)^0, (4.12) 

while 5(1,2) = ^3 = ^(-1,3) = S2 = 5(-2,3) = Si = 0. Then, assuming that the nonreso- 
nance conditions (14.121) are mutually independent, the second-order approximation (14.101) 
becomes 

3 1 ^ 

Xt'' = e J2 + (2^2,(1,-1) + (^2,,E, + ^2,0- ,)E|) 

+ v42,(i,3)EiE3 + A2,(l,_2)EiE_2 + A2,(2,3)E2E3) + C.C. (4.13) 
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By ( I2.17p . the equations for the amphtudes Aij, j = 1, 2, 3, read 



9.A,2 - V^fi(^2)-V,A,2 = -i^^A,3Ai,i, (4.14) 



2 



with = + ^^2, fi(^3) = + ^(^^2) and 

C(i,2) = -4i a2,a sm I y a I sm I y a I sm I y a I - 62- (4-15) 

By (EUl) and (HA2l) we obtain as in Case □ 
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^2,0j) = -P^^lj, j = 1, 2, 3, ^2,(1,-1) = l^ljf, 

^2,(1,3) = ^^Ai,iAi,3, ^2,(1,-2) = ^2,(2,3) = ^^^1,2^1,3 (4.16) 

0(1,3) ^(1-2) 0(2,3) 

with C(p^q) given by (12.191) . while the functions ^2,(1,2), ^2,(1,-3), and ^2,(2,-3) are included 
in ^2,3, ^2,2, and ^2,1, respectively, since (1, 2) = 3, (1, —3) = —2, (2, —3) = — 1 G Af, and 
thus remain undetermined by (12.171) . (I2.18p . 

The system (I4.14p of three nonlinearly coupled transport equations, the three-wave 
interaction equation (cf. |BS90l IKauSOl IKea99l IKRB791 ISW03] ). is the generic system 
describing the macroscopic dynamics of three different modulated pulses, interacting all 
with each other, each being generated by the other two: no subsystem decouples and the 
nonlinearity is symmetric with respect to the amplitudes. 

Ignoring any of the involved pulses by setting, e.g., 3 = 0, leads to the system 

- V^fi(?9l)-VyAi,i =0, 9,Ai,2- V^^fi(7?2)-V^Ai,2 = 0, Ai,iAi,2 = 0. 

This is the case of the non-2-closed system of two pulses which generate via interaction 
a third, discussed in Sec. 14.2.21 Case 1, which, as we saw there, prohibits a complete de- 
scription of the dynamics of the two generating pulses, due to the condition Ai^iAi 2 = 0. 
To obtain the correct dynamics one has to consider (I4.14p . provided of course the system 
of the three pulses is closed. 

Case [3t Only one self-interaction. 

We assume that the resonance condition {2'di,2uji) = {^2,^2) holds, and that no other 
interactions are present. This leads to = 62 = ^(1,-2) = ^-i = 0, and the closedness 
condition is given by the nonresonance conditions 

%i) 7^ 0' j = 2, 3, 5(1,2) 7^ 0, 5(i,±3) ^ 0, 5(2,±3) 7^ 0. 
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The equations (I2.17P yield 



^ a,Ai,3- v^n(793)-v,Ai,3 = 

with C(i^i) given by (14.31) . By (12.181) we obtain the same equations for A2^(jj), j = 2,3, 
^2,(1, 2), ^ 2,(i,±3); ^ = 1,2, v42,(i_i), as in Case [H while v42,(i,i), ^2,(1,-2) are included in 
^2,2, ^2,1 5 respectively, and thus remain undetermined. 

Here, in contrast to (14.141) . the subsystem for the amplitudes 1, Ai^2, given already 
by (14.91) . decouples from the transport equation for 3, since the set {±{'dj,LJj) : j = 1, 2} 
forms itself a 2-closed system, which does not interact with the third pulse. The amplitude 
of the latter travels through the other two unperturbed. Up to relabeling, the full system 
is the closure of the two non-interacting pulses with 5(2,2) = of Sec. 14.2.21 Case 1. 

In this context, we would like to mention the trivial fact that although the resonance 
condition {2-^i,2lui) = (■(?3,a;3) is a special case of {-di +'i?2,'^i + (^2) = {^3,1-^3), where 
{{}i,uji) = (-(92, 1^2)5 the correct dynamics of the first-order amplitudes Ai^i + Ai^2 and ^ 
of the specified approximation X^'^, see (14.10p . cannot be obtained directly from (I4.14p . 
Indeed, in this case the latter gives 



dr{A,,l + Ai,2) - VM^iy^y{Al,l + ^1,2) = -i-^Ai,3(Ai,i + Ai,2) 



drA,-, - V^n(2t9i)-V^v4i,3 = -i 



2A,^A,2 + 



^1,1 + ^1.2 



without the terms in square brackets. However, the latter are included in the cor- 
rect macroscopic equations for the amplitudes Ai^i + Ai^2 and A13 of two pulses with 
(2i?i, 2^71) = (-(93,^73) forming a 2-closed system, cf. (14. 9p . i.e., the subsystem from above. 
The reason for this inconsistency lies in the different form of the underlying multiscale 
ansatz'es (I4.10p vs. (14.51) . Naturally, the correct ansatz for a given (closed) system of 
pulses, is prescribed by the different pulses among them. 

Case m Two self-interactions. 

We assume {2'di,2uji) = (■i?2,i^2) and (2192, 2ti;2) = (Mi^AuJi) = (-(93,^73). Hence, = 
62 = 5(2,2) = ^3 = ^(1,-2) = S-i = '^(2,-3) = 5-2 = 0, and the closedness of the system is 
guaranteed by the nonresonance conditions 

5(3,3) 7^ 0, (5(1,2) = 5(1,-3) 7^ 0, 5(1,3) 7^ 0, 5(2,3) 7^ 0. 

The equations (I2.17P yield 

' drA.^l - VM^y^^yAl,! = -i^^A^A,^2, 



5rA,2- V^n(2t9i)-V,Ai,2 

9,Ai,3- V^n(479i)-V,Ai,3 



. 1 



■(C(1,1)A^ + 2C(2,2)^1,2^1,3), 



^(2^ 2 
' " \^1,2 
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with C(i^i) as in (14.31) and 



C{2,2) = -4i ^ a2,a sin^ i^va) sin (2t9i-a) - 62- 



By ( 12.181) we obtain ^42,(^,3), j = 1,2,3, ^2,(1,-1) as in Case[T], and 

2 

^2,(1,2) = T (C(l,2)^l, 1^1,2 + C(_l,3)^l, 1^1,3)- 

'^(1,2) 



Note that here (1,2) represents also (—1,3), since E1E2 = E3E_i(= Ef). Hence, ^2,(1, -3) 
is included in ^2,(1,2)- The functions ^2,(1,1); ^2,(2,2), ^2,(1,-2), ^2,(2,-3) are included in 
^2,2; ^2,3; ^2,1, ^2,2, respectively, and hence remain undetermined. 

In the present case the third pulse is generated by self-interaction of the second, which 
in turn is generated by self-interaction of the first. This is reflected in the above equations, 
where we observe that (a) Ai^^ = yields ^41^2 = ^1,1 = 0; (b) ^41^2 = implies Ai^i = 0, 
while y4i^3 evolves according to an uncoupled transport equation; and (c) Ai^i = yields 
a coupled system for the dynamics of 2 and 3 of the form (14.91) . cf. the discussion in 
Sec. 14.2.21 Case 2. (In particular, the above system is the closure of the two interacting 
pulses with 5(2,2) = discussed there.) 

Case [5t One self-interaction and one three-wave-interaction. 

We assume (2^91, 2ct;i) = {'1^2,^2) and ('i9i+'i92, co'i+c<j2) = (3'i9i, 3071) = {1^3, u^). This yields 
the resonance conditions = ^2 = 5(i,2) = S3 = '^(1,-2) = = '^(1,-3) = ^-2 = ^(2,-3) = 
5_i = 0, and the closedness condition is equivalent to the nonresonance conditions 

5(2,2) = 5(1,3) 7^ 0, 5(3,3) 7^ 0, 5(2,3) 7^ 0. 

Then, the equations (12.171) yield 

drAi^i - VMi^i)-VyAi^i = -i^^^ {c^A,3A^ + <^A^iAi,2) , 

9.Ai,2 - V^fi(2^0-V,Ai,2 = -i^^ (c(^A,,,A^+^Al,) , 

a.Al,3 - V,fi(3?9i)-V,Ai,3 = -i-^^Ai,iAi,2 



with C(i,i) and C(i,2) given by (14. 3p and ( 14.15^ with '^2 = 2'di and ^^3 = 3^9i, respectively. 
By ( 12.181) we obtain ^2,(3,3), ^2,(2,3), ^2,(1,-1) as in Case[Il and 

^2,(2,2) = T (C(2,2)^12 + 20(1,3)^1,1^1,3). 

0(2,2) 

Here, (2,2) represents also (1,3), since Eg = EiE3(= E^). Hence, ^2,(1,3) is included in 
^2,(2,2)- The functions ^2,(1,1), ^2,(1,2), ^2,(1,-2), ^2,(1,-3), ^2,(2,-3) are included in ^2,2, ^2,3, 
^2,1, ^2,2, ^2,1, respectively, and hence remain undetermined. Note, that the first terms 
on the right hand sides of the amplitude equations are the same as in (I4.14p . However, 
due to the additional self-interaction, we obtain also the corresponding second terms. 
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Here, assuming (a) 3 = 0, yields the system (14.91) under the condition A12A11 = 0, 
cf. Sec. WTTH Case 2 (for 5(1,2) =0). If ^1,2 = we further get A^^i = 0, while if = 
the system reduces to a transport equation for Ai 2- Assuming (b) ^ = 0, we obtain 



- V^^](^9l)■V,Al,l = 0, drAi^s - V,fi(^3)-V,Ai,3 = 0, Al, = ^Ai,iAi,3. 

If 1 = the evolution of 3 can still be described, while ^ = implies 1 = 0. 
Even for ^1,3 7^ these equations yield = Ai^ily^fl^^}^) — V^fi('i9i)) •Vj,Ai,3 on the 
support of v4i 3, which in general implies Ai i = 0. Finally, assuming (c) 1 = 0, we 
obtain the system (14. 7p and 3^1,2 = 0, which in general imply that one amplitude 
vanishes identically allowing to describe the dynamics of the other. Considering the 
succesion in the generation of pulses in the present case, these observations confirm the 
principles postulated in Sec. 14.2.21 

4.3.2 Non-2-closed system of a three- wave-interaction 

We consider a system of three different pulses {±{'dj,Ljj) : j = 1, 2, 3} with uj = > 
0, which satisfies the resonance condition (14. lip , {-di + {^2,^^! + ^2) = {^3,^3), of Case [2] 
in Sec. 14.3.11 but which is not closed under interactions up to order 2. In particular we 
assume that the nonresonance condition 5(2,3) 7^ of (I4.12p is violated, while the others 
are still satisfied. Thus, the pulse {^2 + '^3,^2 + 1^3) = (^^i + 2'd2,oJi + 2u;2) = {'d^^.uji) 
is generated by the interaction of the former two, and the set {^{'dj^ujj) : j = 1, ... ,4} 
consists of four different pulses. 

As we saw in the previous sections, the effect of non-2-closedness on a system of pulses 
results from ignoring the amplitudes of the generated pulses. Thus, assuming that the set 
of the four pulses is 2-closed, we determine next the evolution equations for its first-order 
amplitudes, and discuss subsequently the effects of ignoring the amplitude of ±{{}4,u!/i). 

The second-order approximation (11.130 for z/ = 4 different pulses reads 

4 4 

+ E (^2,o-,^)E,E, + A2,(,,_.)E,E_,) + ^A2,(i,_i)) + c.c. 

l<j<i<4 

with Ej(t, 7) = e^'-'^^*"^'^-' '''^ j = 1, . . . , 4, and when all appearing indices are representants. 
By the assumed interactions we have the resonance conditions 

^(1,2) = ^(1,-3) = ^(2,-3) = (5(2,3) = ^(2,-4) = ^(3,-4) = 0. 

In the case of no further interactions the system of the four pulses is 2-closed if the 
nonresonance conditions 

(5(1,1) 7^ 0, 5(2,2) = '5(-l,4) 7^ 0, 5(3,3) = (5(1,4) 7^ 0, 5(4,4) 7^ 0, 

(5(1,-2) 7^ 0, 5(1,3) 7^ 0, 5(2,4) 7^ 0, 5(3,4) ^ 
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are satisfied. The equations (12.171) yield 







■VyA 






■VyA 






■VyA 


drAl^i 




■WyA 



1 
1 



— 1 



(C(l,2)^l,l^l,3 + C(2,3)^l,3^1,4), 
(C(l,2)^l,l^l,2 + C(^Ai_2Ai^4) , 
-41,2^1,3 



with ^^ = ^1+ ^2, ni'&3) = ^i^l) + ^(^2), ^i = ^2+ ^3, ^{^i) = ^(^2) + ^^(^^3), and 

C(i,2) given by fl4.15p. 



C{2,3) 
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As in the previous sections, the coefficients ^2,(j,i) with (j, i) as in the nonresonance 
conditions can be calculated by the equations (12.181) . while the coefficients ^2,(j,i) with 
(j, z) as in the resonance conditions are included in the second-order amplitudes A2J, 
j = 1, . . . , 4, and thus remain undetermined. 

Note, how the right hand sides of the amplitude equations mirror the two three-wave- 
interactions. In particular, setting identically Ai^^ = 0, the three ffist equations give the 
three- wave-interaction equations (I4.14p . and the fourth equation becomes ^1^2^1,3 = 0. 
As discussed previously (cf. Sec. 14.3.11 Case [2] and Sec. 14.2.21 Case 1), the latter condition 
obstructs the reasonable study of the dynamics of the system (already by restricting the 
choice of initial data) and is safely guaranteed if one of the two amplitudes vanishes. 
However, this leads subsequently to the cases cited above, reducing eventually the system 
for the dynamics of four pulses to a single transport equation for one of them. 



4.3.3 3-closed system of a three-wave-interaction 

We conclude our list of examples by calculating the third-order approximation X^'^ for 
a set of three different pulses {±{{}j,u!j) : j = 1,2,3}, Uj = Q^-dj) > 0, which satisfy 
the resonance condition (14. lip of Case [2] in Sec. 14.3.11 We assume that except for this 
three-wave-interaction, ('(9i-|-^92, t^i+t^2) = {^3,^3), no other interaction (of order 2) takes 
place among the three pulses. Moreover, in order to be able to calculate X^'^ such that 
res(X;f'^) = 0{e'^), we assume that the set is closed under interactions up to order 3. 
According to Definition II. 4[ in the present setting this is the case when the nonresonance 
conditions (or order 2) (I4.12p are satisfied and additionally 

the pairs {k'&i + X'&2, + Aco'2) with 
(«:,A) = (3,0), (0,3), (3,3), (2,-1), (-1,2), (3,1), (1,3), (3,2), (2,3) 
do not characterize pulses, except if either (3'(9i, 3ciJi) = {'1^2, ^^2) or (3'(92, 30^2) = {'di,u!i). 

Hence, excluding the latter exceptional cases (resonances of order 3 within the set of 
pulses), the system is closed under interactions up to order 3 if in addition to (14.120 the 
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nonresonance conditions (of order 3) 

7^ 0' j = 1)2,3, <5(l,l,-2) 7^ 0, 5(2,2-1) 7^ 0, 

5(1,1,3) 7^ 0, 5(2,2,3) 7^ 0, 5(3,3,1) ^ 0, 5(3,3,2) ^ 0. (4.17) 

are satisfied. Assuming that all these nonresonance conditions as well as those in (14.121) 
are mutually independent, the third-order approximation reads 



= Xt'^ + ( ^^3,(1,-1) + E (^3,,E, + A 



2t7I I A T7i2l 



+ ^3,(l,3)ElE3 + y43_(i _2)EiE_2 + ^3,(2,3)E2E3 + ^43^(1^1 _2)EiE^2 + ^3,(2,2,-l)E2E_i 
+ ^3,(l,l,3)EiE3 + y43_(2,2,3)E2E3 + 743_(3_3_i) E3E1 + y43^(3_3^2)E3E2 + C.C.^ (4-18) 

with X^'^ given by (I4.13P in Sec. 14.3. Ij Case [2J There, we determined for X^'^ the 
first-order amplitudes Aij (as solutions of the three- wave-interaction equations (14.141) ) 
and the functions (14.161) . Here, using (14.161) . we obtain from the equations (I2.20p for 
k = 3 the following system of linearly coupled inhomogeneous transport equations for the 
second-order amplitudes A2j: 



drA2,l - V^fi(t9l)-V,A2,l + 1^^(^1,3^2,2 + ^1,2^2,3) 

= + 2^(i--2)IA,2r + 2r/(i,3)|A,3r)A,i - 5,^A,i 

+ ^1,3 7(3,-2)- VyAi;^ + A^2li-2,3yVyAi;3 + ^ 5Z aha^'"^'"" {(^-^ ^Ia) , 



a.^,2 - V^fi(^?2)-V,A2,2 + i^^(Ai,3A^+ ^17^2,3) 

= ^T^^(((r/(2,2) + ^ )|A,2|' + 2r7(2,_l)|Ai,i|2 + 2r7(2,3)|A,3p) A,2 - d'^A^,2 

+ Ai,3 7(3,-1)- VyAi;r + A^-f (-1,3}-^ yAi^s + \Y1 Vj/)^^1,2) , 

aer 

9.^2,3 - VM^3)-VyA2,3 + i^|^ (^1,1^2,2 + ^1,2^2,1) 

= ^T^^(((m3) + ?)IA,3|' + 2r/(3,l)|Ai,ip + 27/(3,2)1^1,2^^,3 - 9^^1,3 

1 . 

+ ^1,1 7(1,2) ■Vj,Ai,2 + Ai,2 7{2,i)-VyAi,i + ai,ae'^'"'{a-VyfAi^. 
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with C(i,2) given by (I4.15p . and 



26^ \cc -)P 

7(p,g) := 2 a2,a{cos{{^p+^q)-a)- cos(?9g-a))a, r^j^p) ■= + h 3cQ-,p,_p). 
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(Note, that the source terms and the coefficients in front of the amphtudes A2J are known.) 
Determining A2J by these equations, and since the set of pulses is 3-closed, we can then 
calculate by fl2.2ip for k = 3 the remaining functions of X^'^, except for the third-order 
amplitudes A-^j. As previously, setting e.g. A^j = 0, we obtain res(X^'^) = 0{e^). 



4.4 Existence of interacting pulses 

We conclude this section of examples by discussing the question whether there really 
exist interacting pulses in lattices. We do this by showing exemplarily the existence of a 
closed system of three different pulses {±{{}j,u!j) : j = 1, 2, 3} which satisfy the resonance 
condition (14. lip (three- wave-interaction) and the nonresonance conditions (14.12^ of Case [2] 
in Sec. l4.3.1l for a one-dimensional lattice (i.e., a chain) of the form (11. 2p with ci = 1, F = Z 
and only a nearest-neighbour interaction potential Vi{x) = y^'^~'<~0{\x\^), where 
in (II. 4p and Va = for a G N \ {1}. The linearized model and the dispersion function 
read, respectively. 



ai(z^+i — 2x^1 + x^_i) — hix^i, 7 G Z, and Vt{'d) = A/2ai(l — cos-^) + 61 > 0, 



where the stability condition (ll.Sp is satisfied iff we assume min{6i, bi + 4ai} > 0. 

For three given pulses ±{i)j,ujj), j = 1,2,3, with ujj := > the resonance 

condition (14. lip is satisfied if and only if the wave numbers '&i,'&2 £ solve the equation 

n\hA+k2A) = {kMA) + k2n{^2))\ ih, h) e (4.19) 

for {ki, /C2) = (1) !)• Moreover, the nonresonance conditions (I4.12p are satisfied if and only 
if (the same) '&i,'&2 do not solve the equations (I4.19P for 

{k,,k2) = {2,0), (0,2), (2,2), (2,1), (1,2), (1,-1). (4.20) 

Inserting the dispersion function ^("(9), equation (I4.19P becomes 

— sgn(ai) cos(A;i'i9i+A;2'i92) = {kl+kl — l)f — sgn(ai)(A;^ cos-i^i + k^ cos^92) 



+ 2A;i/c2 V'(/-sgn(ai) cos'(9i)(/-sgn(ai) cos792) 

with / := ^ + sgn(ai). In particular, for (fci, ^2) = (1, 1), applying the trigonometric 

2|ai| 

1 11-1 1 • • 1 — cos??! , 1 — COS'(92 r/^ n i 

theorem and makmg the substitutions x '■= ^ ; W '■= ^ ^ [u, IJ, the 

equation reads 



sgn(ai)s V'x(l-X)^(l-V') = ^ + sgn(ai)xV' + V (0+sgn(ai)x)(0-Fsgn(ai)^/') 

with s := sgaisiwdi sm'd2) and := ^ > 0. Squaring the left and right hand sides of 

4|ai| 

this equation, we obtain 

-X^(X+V') = — +sgn(ai)0(xV^+X+V^) + (0+2sgn(ai)x^/')V(0+sgn(ai)x)(0+sgn(ai)^/'). 
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Obviously, for oi > there exist no solutions Xji' ^ [0, 1] to this equation. However, for 
fli G (—61/4,0) the stability condition (11.81) still holds, and the equation reads 

502 

9{x,i^) ■= — + x^(x + ^) -(PM + x + i^) + (0 - 2x^) V(0 - x)(0 - V^) = 

with > 1. The function g is symmetric and for 

9ix) ■■= gix, X) = 4x' - 30x' - 30X + 90V4 

we get ^(0) = 9074, ^(1) = 9(0 - 4/3)74, g'{x) = 3(4x^ - 20x - 0), ^"(x) = 6(4x - 0). 
Thus, (7 has a minimum at Xm(0) •= (0/4) (1 + a/1 + 4/0), which lies in (0, 1) if and only 
if < 4/3, and gixM) = -(078) [1 - 12/0+ (1 + 4/0)^/2)] < 0, since for x := 4/0 > 3 
we have (1 + xf^^ > 3a; - 1 <^ (1 + a;)^ > (3a; - 1)^ ^ x{x - 3)^ > 0. 

Hence, for G (1,4/3), i.e. for ai G (-674,-36716) the level set {(x,^) G [0,1]^ : 
5'(X)'0) = 0} is nonempty, and its elements solve 

- s7x(l-x)^(l-V^) = ^ - XV' + V(0 - X)(0 - ^) (4.21) 

with the s = ±1 chosen according to the sign of the right hand side. Changing if necessary 
the sign of "dx or "di in such a way that s = sgn(sin'i9i sin'(?2)! we thus obtain a whole family 
of {^1,^2) G T|, which satisfy the resonance condition fl4.2ip or, equivalently, 

cos(?9i+?92) = / + cos?9i + cos7?2 + 27(/ + cosi9i)(/ + 0081^2) 

with / = 20 — 1 > 1. Among the family of solutions {'&i,-&2) G T| satisfying (14.211) we 
then have to find a pair which satisfies also the nonresonance conditions 

cos(Mi+M2) 7^ ikl+kl-l)f + A;7os79i + cos^^a + 2A;iA;2 7(/+ cos^9i)(/+ cosi^a) 

for the values of {ki,k2) given in (I4.20p . This is equivalent to finding {Xii^) ^ {(X7'0) ^ 
[0, 1]2 : gix^ip) = 0} which satisfy 

X', V-', CVy, XV^, xCCV't^^ with C:=X + V^-0-27(0-x)(0-V^)- 

The intersection of the corresponding curves in [0, 1]2 with the curve of solutions is a 
zero-measure set. Hence, there exist uncountably many {i)i,'&2) G T| which satify (14.191) 
with {ki,k2) = (1,1) (resonance condition) and do not satisfy (14. 19p with the {ki,k2) 
given in (I4.20p (i.e., satisfy the nonresonance conditions). 
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